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Abstract. We present spectroscopic ages, metallicities, and [a/Fe] ratios for 70 globular clusters in M31 that were 
derived from Lick line-index measurements. A new interpolation technique of age-metallicity and a/Fe-diagnostic 
^ . grids is used to account for changes in index strength as a response to abundance-ratio variations, in particular 

' for all of the Balmer-line Lick indices. In addition to a population of old (> 10 Gyr) globular clusters with a wide 

1/^ , range of metallicities, from about —2.0 dex to solar values, we find evidence for a population of intermediate-age 

1/^ ■ globular clusters with ages between ~ 5 and 8 Gyr and a mean metallic ity [Z/H] ~ — .6. W e also confirm the 

presence of young M31 globular clusters that were recently identified bv lBeaslev et all i2004l) . which have ages 
, 1 Gyr and relatively high metallicities around —0.4 dex. The M31 globular cluster system has a clearly super- 

solar mean [a/Fe] — 0.14 ± 0.04 dex. Intermediate-age and young objects show roughly solar abundance ratios. 
We find evidence for an age- [a/Fe] relation in the sense that younger clusters have smaller mean [a/Fe] ratios. 
. From a comparison of indices, mostly sensitive to carbon and/or nitrogen abundance, with SSP model predictions 

CL(' for nitrogen-enhanced stellar populations, we find a dichotomy in nitrogen enhancement between young and old 

M31 globular clusters. The indices of objects older than 5 Gyr are consistent with a factor of three or higher 
in nit rogen enhancement compared to their younger counterparts. Using kinematical data from iMorrison et alJ 
^ ' l|2003) we find that the globular cluster sub-population with halo kinematics is old (Jj 9 Gyr), has a bimodal 

, metallicity distribution, and super-solar [a/Fe]. Disk globular clusters have a wider range of ages, are on average 

more metal-rich, and have a slightly smaller mean [a/Fe] ratio. A cross-correlation of structural parameters 
for M31 globular clusters with spectroscopically derived ages, metallicities, and [a/Fe] ratios shows a correlation 
between half-light/tidal radius and metallicity, which is most likely due to the correlation of half- light/tidal radius 
5^ ' and galactocentric distance. We compare our results for M31 globular clusters with those obtained with the same 

technique for globular clusters in the Milky Way, Large Magellanic Cloud, M81, and other spiral galaxies in the 
Sculptor group. Finally, we compare the globular cluster systems of the two Local Group spirals, M31 and Milky 
Way, with their integrated bulge light. 



Key words. Galaxies: star clusters, abundances, forma- 
tion, evolution. Local Group - Galaxy: globular clusters - 
Galaxies: individual: M31, Milky Way, Large Magellanic 
Cloud, M81, NGC 55, NGC 247, NGC 253, NGC 300 

1. Introduction 

Because they can be resolved into single stars, the stellar 
populations in the Local Group are milestones for under- 
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standing star formation and chemical evolution histories of 
distant galaxies. Globular clusters host simple stellar pop- 
ulations that are characterized by their small dispersion in 
age and chemical composition, and are, compared to the 
mix of stellar populations in galaxies, relatively easy to un- 
derstand. All theoretical predictions for observables of un- 
resolved stellar systems are calibrated on resolved stellar 
populations in Milky Way globular clusters. Since the for- 
mation of globular clusters is temporally linked with ma- 
jor star-formation episodes, globular cluster systems are 
fossil records of the past star forma tion and chemical evo- 
lution histor y of their host galaxies l|Ashman fc ZepilQQSt 
lHarrijl2nnil) . It is therefore of paramount importance to 
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ensure that the predictions drawn from a comparison of 
Galactic and extragalactic globular clusters, and stellar 
populations in general, are not biased by any peculiarity 
of our neighborhood. 

The richest globular cluster systems in the Local 
Group are hosted by the two major spiral galax- 
ies Andromeda (M31; N GC 224) harboring 460 ± 70 
ijBarmbv fc HuchralliooHl and th e Milky Way with 160 ± 
20 globular clusters (|Ha,rriFll2nnil) . Photometry has shown 
that the globular cluster systems of M31 and Milky Way 
have a similar range in metallicity and that both distri- 
butions are bimo dal with peaks around [Fe/H] = —1-4 
and— 0.6 dex fe.g .lBarmbv et al.l2000l:lBarmbv fc Huchral 
1200 ll) . Recently, IPerrett et all l)2002tl spectroscopically 
confirmed this metallicity bimodality. However, M31 hosts 
a larger number of metal-rich clusters, which is reflected 
in the higher m ean metallicity of th e entire globular clus- 
ter system (e.g. iHuchra et al1ll99lh . Although the globu- 
lar cluster lu minosity funct ions are remarkably similar in 
both galaxies (|Ha,rri!jl200lt> . they may diffe r as a function 
ofmetallicity and galactocentric radius ((Barmbv et al.l 
l200l|l . There is also a difference in spatial distribution 
(scale size), where the M31 cluster system appears to be 
more extende d than the Galactic globular cluster system 
(|HaTrijl2nniK Thus, despite the overall similarities of the 
M31 and Milky Way globular cluster systems, there are 
some striking differences in the details. Recently, a revised 
optical/near-IR photometric catalog of 337 confirmed and 
688 candidate globular clusters in M31, that combines 
and homogenizes previous studies, was published by the 
Bologna group and will all ow for more accurate future 
analyses ijGalleti et al.ll20o3. and references therein). 

Integrated-light sp ectroscopy of M31 g lobular clus- 
ters was pioneered by Ivan den Berghl 1)1969,1 who found 
that the M31 globular cluster system extends to much 
higher metallicities than the Milky Way globular clus- 
ter system. Subsequent studies confirmed these results 
Ce.g. lSpinrad fc Schweizerlll972l:lRabinlll98l[l . A thorough 
study, based on spectroscopic Lick line indices of 19 M31 
and 17 Galactic globular clusters, bv lBurstein et al]l)l984(l 
confirmed this and showed that M31 clusters have sig- 
nificantly stronger H/? and CN absorption indices at a 
given metallicity. After carefully checking the influence of 
hot horizontal-branch stars and/or blue stragglers, and 
younger ages on these indices, the authors concluded that 
age is the most favorable, but not sufficient, explanation 
for the systematic differences. They speculated that their 
sample was dominated by relatively young disk globular 
clusters, which could explain the offset in H/3, but not the 
strong CN indices. The possibility of an enhanced car- 
bon and / qr nitr ogen abundance was considered unlikely. 
iBursteii] l|l987t) later showed that K giants contribut- 
ing to the integrated light of M31 clusters must have at 
least as strong CN absorption bands as the most extreme 
Galactic giant stars. In a subsequent analysis of the Call 
H-l-K absor ption fea t ure o f nine metal-rich M31 globu- 
lar clusters, iTr iniccnl ^^E^ rejected the hypothesis that 
hot horizontal-branch stars and/or blue stragglers boost 



the H/? index of the integrated light. Modeling the spec- 
tra of M31 clusters he also found that the strong CN in- 
dex is consistent with roughly a factor of 10 in CN excess 
in both dwarf and giant stars. iDavidgel l|l990al) obtained 
near-infrared spectra of four luminous M31 globular clus- 
ters and compared them with Galactic counterparts. He 
found a clear enhancement in the 2 /xm CN band and 
concluded that the increased CN abundance in M31 glob- 
ular clusters is likely to be of primordial origin. From a 
comparison of CO absorption bands in M31 clusters and 
the integrated light of M32, whi ch is known to host an 
intermediate-age population fe.g. lDavidgjll990b|) . he rea- 
soned that, if both stellar populations have similar metal- 
licities, luminous M31 globular clusters may contain an 
intermediate-age component. 

In a comprehensive spectroscopic investigation of 149 
globu lar clusters in M31 , using specifically defined line in- 
dices, iBrodie fc Huchral l)l990(l confirmed a CN offset rel- 
ative to Milky Way globular clusters, but failed to iden- 
tify a stronger mean H/3 index, although good agreement 
was found between the measurements of matched objects. 
This fact suggests that sample selection in the M31 clus- 
ter system might be playing an important role and implies 
the presence of multiple globular cluster sub-populations 
with different mean H /? indices and pote ntially different 
chemical compositions. iBrodie fc Huchral reported also an 
enhancement in their Can H-t-K index for M31 globular 
clusters, at a given metallicity. Assuming that the CN en- 
hancement is due to an increased nitrogen abundance they 
identify aluminum satellite lines as responsible for an in- 
creased H-l-K index, because of the correlation of nitro- 
gen with aluminum abundan ce, at least in metal-rich stars 
l)Norris fc Pilachowskil ll985'). 

An alternative explanation for the observed CN off- 
set is lower mean su rface gravities. a. s these elevate both 
CN and Can H-l-K l|0'Connel]l 11971 . The surface grav- 
ity can be lowered by decreasing the age; however, other 
indices such as Mg2 (a measure of the Mg6 triplet and 
the magnesium hydride band, MgH) increase significantly 
as age is decreased, and such an increase is not observed. 
However, the Mg2 index can be kept constant if the in- 
trinsic magnesium abundance and age are both decreased 
appropriately. The former could be a result of a low 
[a/Fe] ratio at constant total metallicity. This is expected 
for stellar populations which formed late fror n material 
enriched by type II and la su pernovae (e.g. iMatteuccl 
11994 iMatteucci fc Recchil EoOlh . Later work using near- 
UV spectra analyzed the strength of the NH band at 3360 
A (closely linked to th e nitrogen abundance) for several 
M31^ globular clusters llPonder et alJll998l: Ei fc Bursteii] 
'2003': Bu rstein et alJl2004f) . These^udies confirmed that 
M31 globular clusters are indeed en hanced in nitrogen. 

Furthermore, previous studies llBarmbv et al.l l200(]l: 
iBeaslev et allEoMlBur stein et al ■l2004l) showed the pres- 
ence of young M31 globular clusters. These young objects 
have spectra which resemble those of young globular clus- 
ters in the Large Magellanic Cloud with ages around 1 
Gyr. Some of these M31 clusters appear to belong to a 
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Table 1. Radial coverage of background sectors. r„iin, 
'"max, and r,„odian are the radii of the fiber pointings with 
the minimal and maximal galactocentric distance, and the 
median distance of the entire sub-sample. 



sector 


fmin [kpc] 


rmax [kpc] 


r-median [kpc] 


zO 


0.55 


0.85 


0.60 


zl 


1.56 


4.96 


2.48 



distinct su b-population of t hin d isk globular clusters dis- 
covered by iMorrison et al.l l(2004l) . 

It is obvious that a study of ages and chemical compo- 
sitions for a large sample of M31 globular clusters is ur- 
gently needed to resolve some of the above issues. Previous 
studies have shown indications of young ages for some 
M31 globular clusters and clear signatures of chemical dif- 
ferences between the globular cluster system in M31 and 
other galaxies. 

In this paper, we focus on the analysis of spectro- 
scopic features of M31 globular clusters using Lick line 
indices from which we derive ages, metallicities, and 
abundance ratios. Newly available simple stellar popula- 
tion (SSP) models for wel l-defined chemical compositions 
(IThomas et al ■120031120041) make these parameters accessi- 
ble with unprecedented systematic accuracy. We use high- 
quality Lick index measurements and a new iterative pre- 
scription to derive these parameters that maximally re- 
duces systematic uncertainties in the internal calibration 
of Lick indices. 

The paper is organized as follows: ^ describes the 
dataset and its calibration to the Lick system. In the fol- 
lowing Section fj^lwe describe our method and perform a 
consistency check for ages and metallicities derived from 
different diagnostic grids. Results on M31 globular cluster 
ages, metallicities, and abundance ratios are presented in 
Section ^ along with a comparison with globular cluster 
systems in other galaxies. Our findings are discussed in ^ 
and summarized in fJHl 

2. Data 

2.1. Globular Cluster Spectra 

We use spectroscopic observations obtained with the Wide 
Field Fibre Optic Spectrograph fWYFFOS. Bridges..l998i) 
at the 4m William Herschel Telescope on La Palma, 
Canary Islands. Th e dataset is described in detail by 
IPerrett et al.l l|2002fl and we refer the reader to this pa- 
per for a discussion of technical details. Figure ^ shows 
the field of view with all fiber positions. Open symbols in- 
dicate globular clusters; crosses mark the position of back- 
ground fibers. Sky corrections for globular cluster spectra 
were performed by subtracting a mean sky spectrum that 
was constructed from several nearby fibers placed in close 
vicinity to each individual globular cluster. Throughout 
the following analysis we use these sky-subtracted spec- 
tra. 



2.2. Spectra of the Diffuse Light 

In addition to our globular cluster measurements (see cir- 
cles in Fig. nj, we also obtained 207 integrated spectra 
of the diffuse stellar population at different galactocentric 
radii. Each individual spectrum has too low a signal-to- 
noise (S/N) to allow reliable Lick index measurements. In 
order to maximize the S/N ratio, we select spectra from 
the inner Bulge field (within 5 kpc projected galactocen- 
tric distance) and subdivide this background sub-sample 
into an inner and outer sector, zO and zl, respectively. The 
radial sampling of the combined two sectors is described 
in Table Individual spectra in each sector are then co- 
added. Before summing the spectra, a careful smoothing 
to the Lick resolution of each single spectrum was per- 
formed, accounting for the varying velocity dispersion at 
each individual fiber position and the changing spectral 
resolution as a function of chip position. Velocity disper- 
sions as a function of galactocentric dist ance were cal- 
culate d with the disk-bulge-halo models of IWidrow et alJ 
l(2003l) . The average S/N ratio of both combined back- 
ground spectra exceeds ~ 25 per A. 

2.3. Calibration of Lick Indices 

Prior to the measurement of Lick line indices, our spectra 
were broadened to the characteristic Lick resolution. The 
steps undertaken du ring this pr ocedure were adopted from 
Worthev k Ottavi ani ( 199(1) : iKuntschneil fl OgS. 200dh; 
Vazdekis fc Arimota (,1999il : IPuzia et all l|2002bi. i2004ali . 
and IPuzial l|2003bl) . A particular characteristic of the 
WYFFOS instrument is its resolution as a function of 
chip position. To account for this, we computed 2-dim. 
resolution maps using arc exposures that were taken close 
in time with the science frames. The resolution maps 
were then used to smooth o ur spectra to the differential 
Lick resolution of ^ 8 A (see lWorthev fc Ottavianilll997t 
for details). We tested this transformation by matching 
our globul ar cluster list with th e M31 globula r cluster 
sampl e of IWorthev et al.l l|l994|) and iKuntsch ner et alJ 
(Hooi). The comparison gives six matches (B012-G64, 
B015, B193-G244, B218-G272, B225-G280, B311-G33) 
that cover a range in [Fe/H] from —2.0 to —0.3 dex. Small 
linear correction terms were found for all indices, that in 
all cases had a significance of ^ 1 cr (see Table . These 
systematic offsets were applied to all subsequent Lick in- 
dex measurements^. 

Blue Lick indices (from CNi to Fe4383) were measured 
on spectra taken with the blue arm of WYFFOS, which 
covers the spectral range ~ 3700 — 4500 A. The remain- 
ing Lick indices (Ca4455 to Fe5709) were measured on 
spectra taken with the red-arm with a spectral coverage 
from - 4400 to 5800 A. Spectra with S/N of - 25 per A 

^ Note that bad pixel masking leaves two out of the six spec- 
tra useful for calibration of the H7A index. The two remaining 
good-quality spectra are in very good agreement with the Lick 
standard measurements. However, we warn that the rms of this 
calibration cannot be trusted given the small sample. 
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Fig. 1. Positions on the sky of studied globular clusters (open symbols) and background fields (crosses). Globular 
clusters are parameterized by their age and metallicity. Squares indicate globular clusters with ages below 5 Gyr, 
circles show intermediate-age globular clusters with ages between 5 and 8 Gyr, while diamonds mark the position of 
old globular clusters {t > 8 Gyr). Old objects populate the entire studied field. Intermediate-age globular clusters seem 
to prefer the region within the central ~ 5 kpc, while young clusters are preferentially located in the outer disk. The 
metallicity is indicated by color-coded symbols. Red symbols show metal-rich globular clusters with [Z/H] > —0.8, while 
blue symbols show their metal-poor counterparts. Note the clear concentration of metal-rich objects around the nucleus 
of M31. Background spectra were grouped in sectors zO and zl according to their galactocentric distance indicated 
by the two inner circles. A few globular clusters (e.g. MII=G1) which are located at large galactocentric distances 
fall off the image, similar to some background pointings. The underlying im age was taken from DSS2/red. Shaded 
trails indicate stella r streams discovered by the INT/CFHT collaboration fsee llbata et alJli'oOlUFerguson et alJl2002l: 
iMcConnachie et alJEo0 41. While the southern stream is the most massive (~ 10* — 10^ M©), the northern spur carries 
abou t an order of magnitude less mass and appears to be kinematically associated with NGC 205 ( McConnachie et all 
Reticles mark pointings for which deep ph otometry is availab le. Going f rom sma ll to larger galactoce ntric 
distances the fields are labeled as S OI (..Saraicdini fc Van Duvnell200l|) . B03 (Brown et all [2003,) . DOl (^urrcUeiai] 
1200 ll) . and D04 llDurrell et alJl2004l) . 
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Table 2. Correction terms of the transformation to 
match the Lick/IDS instrumental resolution. The correc- 
tions were used to tr ansform our index me asurements to 
the Lick/IDS system IWorthev et alJ l)l994j) . in the sense 

-^Lick -^measured ~t~ ^- 



Index 




c 


rms 


TT r 

HdA 


— 


127 


2 


041 




—0 


437 





697 







001 





096 







051 





084 


Oa4zz7 


— 


407 





956 




— u 


Q 1 K 

3to 


3 


059 


H7A 


— 


488 





514 


H7F 


— 


094 


1 


106 


TT1 /I 

]<e4383 





764 


2 


480 


Ca4455 





541 





796 


Fe4531 


1 


086 


1 


445 


Fe4668 





007 


2 


193 


H/3 





148 





175 


Fe5015 


-0 


300 


1 


547 


Mgi 





004 





025 


Mg2 





025 





041 


Mgfe 


-0 


123 





695 


Fe5270 





730 





589 


Fe5335 





243 





576 


Fe5406 


-0 


064 





549 


Fe5709 





355 





489 



or less were rejected from the subsequent analysis. As we 
propagate the full error budget of the spectra into the line 
index measurement, some indices show unusually large er- 
ror bars due to masked bad pixels and/or cosmics. These 
contaminated index measurements were also rejected from 
the following analysis steps. With these data quality re- 
quirements 29 globular clusters remain in our sample. 

2.4. Data from other Sources 

To maximize the statistical significance of the glob- 
ular cluster sample, we includ e five M31 g l obular 
clusters from the Lick group ijWorthev et alJ 
IWorthev fc Ottaviani"l997^ and 36 globular clusters from 
the |j3easley-et al, (2004) sample. We note that all index 
measurements are performed in the same passband sys- 
tem as the SSP models, in particular the fitting functions 
that are used to compute theoretical predictions. 

We find two matches between our globular cluster data 
and the Beasley et al. dataset. The Lick index measure- 
ments for the two clusters, i.e. B225 and B243, are in very 
good agreement (within 1-cr), especially for the widely- 
defined molecular-band indices (e.g. CN2 and Mg2, but 
not for G4300^). Very good agreement is also found for 
indices crucial to further analysis steps, such as Balmer, 



iron, and magnesium indices: US a, Hja, H/3, Mg2, Mg6, 
Fe5270, Fe5335, etc. Because of the smaller statistical er- 
rors of the Beasley et al. data we adopt the latter for the 
following analysis. Hence our final dataset contains high- 
quality spectra (S/N ^ 25 A~^) for 70 globular clusters 
in M31. The full set of index measurements is presented 
in Tables IXTl and IX!2l in the Appendix. Not e that object 
B166 is suspected to be a foreground star ijCalleti et alJ 
However, since this is not certain, and because we 
do not find any features in our spectrum that would be 
characteristic of a foreground star, we mark this object as 
potentially misclassified but keep it in our globular cluster 
sample. Note that based on their colors and radial veloc- 
ities, all other objects were classified as gen uine globular 
cluste rs by several independent studies (see iGalleti et alJ 
I2OO4I) . Figure |21 shows representative spectra from our 
sample. 

We do not include indices from iBurstein et all 
since they were measured in the old Lick passband sys- 
tem, for which no model predictions with variable abun- 
dance ratios are available. Howe ver, we find eigh t matched 
M31 clusters in our and the Trager et alJ l)l998l) data, who 
re-measured Lick ind ices for M31 an d Milky Way globu- 
lar clusters from the iBurstein et al.l sample with slightly 
modified passband definitions compared to the ones used 
in this work^. We therefore briefly investigate the pre- 
viously detected H/3 of fset between Galactic and M31 
globular clusters in the IBurstein et al.1 l)l984(l study. In 
Figur e Bl we reproduce the H/3 and CN2 vs. Mg2 plots 
from IBurstein et al.1 l)l984|) . We include oi ily the eight 
M31 globular clusters that overlap with the iTrager et alJ 
data and were re-observed in this study (soH d circles) and 
augme nt the sample with clusters from the iTrager et alJ 
l|l998|) study that meet our quality selection criteria (five 
open circ les). Also included a re 12 Galactic globular clus- 
ters fromlPuzi^^g^^lJ 1 2QQ2dl (solid squares), and 13 clus- 
ters from ITrager et al.l (jl998r (open squares). 

We find four match es b etween the G alactic cluster 
samples of IPuzia et"al] and ITrager et al.1 and plot for 
matched objects the indices measured bv lPuzia et al 

1 be- 
cause of their smaller error s. The systematic offsets for the 
four matches in the sense IPuzia et al.l - ITrager et al.l are 
AH/3 = 0.31 ± 0.28 A and ACN2 = -0.02 ± 0.04 mag. 
The matched objects are marked with large open squares. 
From the eight matc hed M31 glob ular clusters between our 
sample and that of ITrager et al.l we find in general very 
good agreement between the two index measurements for 
all relevant indices, such as Balmer, Mg, and Fe indices. 
In particular, the offsets, in the same sense as above, are 



^ We find a ~ 1.5 A systematic offset towards higher values 
in the Beasley et al. sample, although one of our measurements 
has a very large error of 1.2 A. 



^ The difference in passband definitions b e tween 
the ones adopted in th is work (IWorthev et_alj _|l99^ 
IWor thev fc Ottavianil Il997^ and the ones of iXraeer e^j] 
■ 19981*) is significant for some indices. Linear tra nsformations 
between the two systems were computed by IPuzia et alJ 
i2004al) . H/3 has identical definitions in both systems. For 
the CN2 the offset is ACN2(Tr98 - Wo94)0.014 ± 0.020 mag. 
Where applicable, these linear corrections are applied to the 
Trager et al. data. 
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Fig. 2. Representative spectra from our M31 globular cluster sample for selected targets, sorted by metallicity (in- 
creasing from top to bottom). Blue and red spectra are plotted separately. Derived ages (in Gyr), metallicities, and 
[a/Fe] ratios from the subsequent analysis steps are indicated. 



AH/3 = 0.06 ± 0.13 A and ACN2 = 0.026 ± 0.023 mag. 
The systematic offset in H/3 is clearly less significant than 
the one found between the Galactic cluster samples. 

Although the error bars of the M31 sample are rela- 
tively large, there is no significant systematic offset be- 
tween M31 and Milky Way globular clusters in the newly 
acquired data (solid symbols). Fitting straight lines to the 
new data in the H/3 vs. Mg2 plot yields the following linear 
relations: 

M31 : H/3 = (-4.05 ± 1.03) • Mg2 + (2.58 ± 0.18) 
MW : H/3 ^ (-3.95 ± 1.50) • Mg^ + (2.78 ± 0.22) 

The rms for both relations are 0.30 and 0.21 A for M31 
and Galactic globular clusters. Since the gradients are vir- 
tually identical, the zero-point offset is a measure of the 
systematic offset. It is smaller than the mean rms of both 
relations. If we split the datasets at Mg2 = 0.2 A, the 
systematic offsets for both high and low-Mg2 subsets dis- 
appear, i.e. they are less significant th an ^ 0.5g. Ho wever, 
if the Galactic globular clusters from iTrager et alJ (small 
and large open squares) are considered, we confirm a 
systematic offset AH/3 « 0.5 A (^ 3a) at Mg2 > 0.1 
mag between M31 and Milky Way (see also Fig. 5k in 
iBurstein et"al]ll984 . In summary, a combination of sam- 
ple selection, small number statistics, and a potential cal- 
ibration inconsistency might have introduced spurious H/3 



offsets between Galactic and M31 globular clusters in the 
IBurstein et al. work. 



The new data do not appear to alter the conclusion 
of IBurstein et al.] on the CN offset. However, for metal- 
rich globular clusters the offset in CN seems smaller than 
for metal-poor globular clusters. In general, this exercise 
demonstrates that for drawing reliable conclusions from a 
direct comparison of datasets, a homogeneous and consis- 
tent sample selection of these datasets is crucial. 



Later in th is work we show that the speculation of 
IBurstein et al.l that the H/3 offset was produced by a sys- 
tematically younger mean age of their observed M31 glob- 
ular cluster sample is likely to be the result of Milky Way 
globular cluster sample selection. From their sample of 19 
M31 globular clusters, 16 had their indices on high-quality 
spectra re-measured in subsequent studies. In fact, in this 
study we find that the mean age of this M31 sub-sample is 
^10 Gyr with only four globular clusters between 6 and 9 
Gyr. However, their intuition that M31 might host young 
globular clusters turns out to be correct, if a larger more 
representative sample of M31 globular clusters is consid- 
ered. 
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Fig. 4. Comparison of diagnostic grids for Lick in dex measurem ents, parameterized for different [a/Fe] abundance 
ratios and ages. Model predictions were taken from ixhomas et ak (.2003, 2004) . The left and middle panel show two 
grids constructed from [MgFe]' as a metallicity indicator and two different Balmer line indices as age indicators, i.e. 
H7A and H/3, respectively. The plots show models grids for two different values of [a/Fe] — 0.0 (solid lines) and 0.5 dex 
(dashed lines). Note the increasing differential offset between the grids towards higher metallicities. The right panel 
shows the [a/Fe] diagnostic diagram, which features the Mg2 index as a representative a-element indicator and <Fe> 
as a iron-abundance proxy. The diagram shows the dependence of iso-[a/Fe] as a function of age. Dotted lines show 
the model predictions for 3 Gyr old populations and solid lines for a 13 Gyr counterpart. 



3. Analysis 

3.1. Comparison with SSP Models 

In the following Section we compare line index measure- 
ments for M31 globular clusters and background fields 
with predictions of SSP models, which were calculated for 
stell ar population s with well-defined [a/Fe] abundance ra- 
ti os llThomas et al 2003 . 2Q p4). These models are based 
on lWorthev et all |l99l and lWorthev fc Ottavia^lll997tl 
fitting functions, response functions from ^"^ ^^j^col^Bell 
and Korn et al. (2004, in preparation), solar-scaled 



stellar evolutionary tracks, an d employ the princijples o f 
the fuel consumption theorem ijRenzini &: B uzzoni"l986'). 

Spectroscopic ages and metallicities of unresolved glob- 
ular clusters are traditionally derived from Balmer line 
indices and a com bination of m etal indices. According to 
FPuzia. (,2003atl and lPuzia et al.i (2004b 'I the best combina- 
tion of Lick indices to construct an age-metallicity diag- 
nostic diagram depends on internal uncertainties of the 
Lick system and the data quality. In our case the most 
reliable combination is II7A and [MgFe]'^. The combina- 
tion of magnesium and iron-sensitive indices, in the form 
of the composite index [Mg Fe]', is an excellent metallicity 
tracer l)Puzia et alJl2002m and shows the weakest sensi- 
tivity to [a/Fe] abundance ratio variations ^Thomas et al.l 
l2003tl . Among all Balmer Lick indices, H7A, followed by 
H/3 and H(5a, are least affected by internal systematics, 
such as systematic uncertainties in the fitting functions. 
However, we note that all Balmer indices are influenced by 
abundance ratio variations (e.g. a/Fe), in particular the 
higher-order Balmer indices. Although response functions 
exist only for a selected set of elements (Trioicco fc_Bcll 



[MgFe]' = ^Mg6 (0.72 • Fe5270 + 0.28 • Fe5335). 



[l993), with the new SSP model predictions these varia- 
tions are under control to first order. 

Considering the signal-to-noise of a dataset, the accu- 
racy of the calibration to the Lick system, and internal sys- 
tematic uncertainties of Balmer Lick indices, IPuzia et alJ 
l)2004al) computes a reliability parameter TZ, that indi- 
cates the quality of a Balmer Lick index as an age indi- 
cator. Based on this classification we have decided to de- 
rive ages and metallicities from diagnostic diagrams con- 
structed from the Balmer indices H7^, H/3, and H^^i, only. 
The narrowly-defined Balmer indices, H7i? and ^5f, have 
low signal-to-noise and measure the same absorption fea- 
tures as H7/1 and ^5 a (hence, do not provide additional 
information). They have the lowe st reliability param eter 
TZ among all Balmer indices (see iPuzia et all 120044 - so 
that we refrain from including them in subsequent analy- 
ses. 



3.2. Fitting Routine 

To derive reliable ages and metallicities, we need to utilize 
a diagnostic grid from which the [a/Fe] ratio of the ob- 
served stellar population can be derived. We tested several 
combinations of Lick indices and, based on the internal 
systematics, the relative separation of iso-[a/Fe] tracks, 
and the mean measurement error, decided on an [a/Fe] 
diagnostic grid that is constructed from the indices Mg2 
and <Fe>. The Mg2 index measures the strength of the 
Mg6 and MgH features at - 5100 A. The mean iron index 
<Fe> is a composite index built from the indices Fe5270 
and Fe5335, i.e. <Fe>= VFe5270 -I- Fe5335, which trace 
two strong iron lines at similar wavelengths (A 5270, 5335 
A) as Mg2. However, since the [a/Fe] diagnostic grid it- 
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Fig. 3. H/3 and CN2 vs. Mg2 plots for M31 and Galactic 
globular clusters. Milky Way globular clusters are indi- 
cated as squares, while M31 globular clusters are marked 
by circl es. Solid squares show data taken from lPuzia et alj 
ll2nn2bt) Open squ ares a re Galactic globular clusters 
from iTraeer et all l)l998() . Large ope n s quares illus- 
trate matches between the IPuzia et al.l and iTraeer et al.l 
datase t. Solid circles are M31 globular clusters included 
in the iBurstein et alJ dataset that were re -measured in 
this st udy. Open circles add M31 data from ITraeer et alJ 
lll998l). See text for details. 



self is degenerate in age and metallicity, the best possible 
determination of ages, metallicities, and [a/Fe] ratios re- 
quires an iterative fitting of the data, feeding information 
acquired from the a/Fe-grid into the age/metallicity-grid. 

Extensive tests have shown that the values derived 
from both diagnostic diagrams converge after a few it- 
erations. Our routine starts with the determination of 
[a/Fe] ratios from the Mg2 vs. <Fe> diagram, assuming 



a 13 Gyr old population, and is based on an minimiza- 
tion technique. The age and metallicity of each globular 
cluster is then determined with the same minimiza- 
tion approach on an individually interpolated^ H7A, H/3, 
and H(5a vs. [MgFe]' grid for a given [a/Fe]. For glob- 
ular clusters with index measurements that do not fall 
on the diagnostic grid, we assign the most likely value, 
extrapolated along the most likely error vector which is 
constructed from the square sum of individual index un- 
certainties for each globular cluster in a given diagnostic 
grid. The age and metallicity values are then fed back 
into the [a/Fe] diagram and the process is iterated. In 
most cases this process requires only a few cycles for 
the values to converge. This method is superior to any 
age-metallicity determination based on the age-metallicity 
grid only. As shown in FigureQ] the isochrones of the age- 
metallicity grid can change significantly up to a few typical 
1-a measurernent errors a s a function of [a/Fe] (see also 
iThomas et al.l2nn3[EnO^ . To control these effects and ob- 
tain meaningful age/metallicity results it is imperative to 
include abundance ratio variations in each age-metallicity 
determination based on Lick index measurements. 

Ages, metallicities, and [a/Fe] ratios for all M31 sam- 
ple globular clusters are summarized in Table IA.3I in the 
Appendix. 



3.3. Consistency Check 

In FigureElwe compare age and metallicity determinations 
from the three diagnostic grids using different Balmer- 
line indices, ^Ija, Hp, and 116a, as age indicators (see 
Fig. IHJi. The rms scatter in the comparison of ages de- 
rived from two different grids varies between 3.6 and 4.7 
Gyr. Note that the age scatter increases from young to 
old ages, which is an effect of an increasing age resolu- 
tion with decreasing absolute age. Given the quality of 
our dataset, the rule of thumb for the average age uncer- 
tainty is At/t K 1/3. The consistency check shows that the 
rms-scatter in all upper panels of Figure|Slis in most cases 
comparable or larger than the mean measurement error. 
This is an indication for the presence of systematic offsets 
in the absolute age scale between two diagnostic grids. 
The rms-scatter in all metallicity consistency plots agrees 
very well with the mean measurement error, and shows 
the good absolute metallicity calibration of individual di- 
agnostic grids. It is not too surprising that the metallicity 
values derived from different grids agree very well for in- 
dividual globular clusters, since all plots use [MgFe]' as 
metallicity proxy. A good estimate for the mean metallic- 
ity uncertainty is 0.2 dex. 

The consistency shown between the model predictions 
should be taken as a demonstration of the quality of 
age/metallicity predictions that arc feasible with today's 
high-quality datasets and SSP models. Rather than being 
an endorsement of the particular model used in this study, 



^ We linearly interpolate the grid in age and Z-space 
metallicity log-space) . 



not 
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Fig. 5. Comparison of age and metallicity predictions from the three different age-metalhcity diagnostic grids shown 
in Figure El Upper panels show the comparison of ages; lower panels illustrate the comparison of metallicity values. 
The rms scatter is indicated in the upper left corner of each plot and is the scatter around the one-to-one relation. 
(A/) indicates the median error for this panel. In the upper panels median error bars are indicated at the abscissa 
and ordinate, calculated as the sample median of the averaged globular cluster asymmetric error. Open circles mark 
outliers which were not included in the rms determination. 



we point out that SSP model predictions for stellar popu- 
lations with well-defined abundance ratios are essential to 
reach the same quality of age/metallicity determinations 
as presented here. Other SSP models which do not take 
the systematics of abundance-ratio variations into account 
may lead to similarly good fits for a given diagnostic grid, 
but the results will suffer from systematic uncertainties 
described earlier. However, even for state-of-the-art SSP 
models the litany of systematic uncertainties is long^. 

In general, we find good congruence between age 
and metallicity predictions between diagnostic plots us- 
ing the three indices H^^, H^a, and II/3. The relative 
age accuracy (~ 3 — 5 Gyr) allows us to disentangle old 
from intermediate-age and young globular cluster sub- 
populations within a wide range of metallicities. 

^ Systematic uncertainties are still hard-wired into all cur- 
rent SSP Lick-index predictions and might, if accounted for, 
significantly improve future age, metallicity, and abundance 
ratio estimates. To name a few "to-dos": control of system- 
atic uncertainties introduced by variations in horizontal branch 
morphology, self-consistent use of enhanced stellar evolution- 
ary tracks/isochrones in SSP models, extension of response 
functions to more elements, self-consistent calculation of in- 
ternal energy production for various chemical compositions, 
etc. 



4. Results 

4.1. Ages and Metallicities 

FigureElshows the M31 globular cluster data, indicated by 
solid circles, toge ther with theo r etical SSP m odel predic- 
tions taken from iThomas et alJ l)2003l l2004l) . The model 
grid covers ages from 1 to 15 Gyr and metallicities [Z/H] 
from —2.25 to -1-0.67 dex, indicated by sohd and dashed 
lines respectively. The 5 Gyr isochrone is indicated as a 
thicker line to guide the eye when comparing diagnostic 
plots constructed with different Balmer-line indices. We 
note that direct comparisons of data and diagnostic grids 
are misleading when the [a/Fe] abundance of each individ- 
ual globular cluster is not taken into account (see previous 
Section). However, to allow an assessment of the mean age 
and metallicity distribution of all globular clusters and to 
check for consistency of these predictions from all three 
diagnostic diagrams, we plot the data together with di- 
agnostic grids for stellar populations with [a/Fc] = 0.3 
dex. 

The majority of our sample of M31 globular clusters 
are consistent with old ages. However, the sample spans 
a wide range of ages, with ~ 10% of the entire sample 
younger than 5 Gyr. Several globular clusters have index 
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Fig. 6. Diagnostic Plots: Hja, H/3, and USa vs. [MgFe]' age-metallicity diagnostic plot for stellar populations with 
[a/Fe] = 0.3 dex. The lower right panel shows the fa/Fe] diagnostic plot with Mg2 vs. (Fe). Solid circles are M31 globular 
clusters. SSP model predictions were taken from Thomas et al. (2003, 2004) and are shown in the age-metallicity grids 
as solid isochrones for stellar populations with ages 1, 2, 3, 5, 10, and 15 Gyr and as dashed iso-metallicity tracks for 
metallicities [Z/H] = —2.25, —1.35, —0.33, 0.0, +0.35, and +0.67 dex, for a stellar population with [a/Fe] = 0.3 dex. 
The SSP model in the lower right panel shows iso-[Q;/Fe] tracks for metallicities [Z/H] — —2.25 to +0.67 and an age of 
13 Gyr. Open circles mark M31 bulge fields at different galactocentric radii, see Tabled NB87, which has the highest 
metallicity among our sample globular clusters, is the only object with [MgFe]' « 4 A and (Fe) « 3 A. It should be 
stressed again that the shown diagnostic grids are not suited to compare individual globular clusters, because of the 
a/Fe-sensitivity of Balmer indices. All grids have to be interpolated for each individual globular cluster according to 
its [a/Fe]. However, H/? shows the weakest response to a/Fe variations (see Fig. 0J. 



values which are consistent with formal ages similar or less 
than ~ 1 Gyr at an average metallicity ([Z/H]) ft! —0.4 
dex. Unfortunately, the models do not cover ages below 1 
Gyr and we do not attempt to extrapolate the grid, be- 
ing aware of the potential uncertainties introduced by the 
AGB-phase transition which occurs at ages at or below 



~ 1 Gyr. Depending on the total metallicity, the AGB 
can contribute ~ 50% to the total bolometric light at 
these young ages ijRenzini fc Buzzonilll986t iFerraro et alJ 
Theoretical models show that, although the frac- 
tional contribution of the AGB to the total light is sub- 
ject to only small model-to-model fluctuations, the exact 
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onset and duration of the AGB-phase transition depends 
on the treatment o f classical mixing vs. overshooting (e.g. 
lAlongi et al.1ll99.ll) . 

Cluster s with such young ag es are consistent with the 
findings of iBeaslev et al.l ((2004 ) who derive ages for their 
sub-sample of young globular clusters from a comparison 
with empirical spectra of young LMC globular clusters. 
This comparison is backed by their H/3 measurements be- 
ing consis tent with ages betw e en 10 and 800 Myr within 
the latest iBruzual fc Charlol l)2003 ) SSP models^. 

We find several M31 globular clusters with near-solar 
metallicities and one object (NB87) reaching a formal 
[Z/H] « +0.5 ± 0.1 dex at a formally derived age of 9 ± 2 
Gyr, and a [a/Fe] — 0.09 ± 0.05. It is the most metal-rich 
globular cluster in our sample. Because of its relatively 
high metallicity, this object might host a stellar popula- 
tion representative of those in massive early-type galaxies, 
although it is slightly under- abundant in a-elements. 

4.1.1. Age-Metallicity Distributions and Correlations 

In the following we compute mean ages and metallicities 
for our sample globular clusters using the information ac- 
quired from all three diagnostic grids. Based on the 7?.- 
parameter ranking and th e sensitivity of Balm er lines to 
abundance variations (see iPuzia et al.ll2004a|) . we assign 
a weighting of 1.0 to values obtained from the Il7^ grid 
and a weighting factor of 0.5 to results obtained from the 
H/3 and RSa grids. The reason for this weighting pattern 
are internal systematics in the age-metallicity diagnostic 
grids utilizing diffe rent Balmer-line indices (for details see 
IPuzia et al]l2004a|) . 

Using this averaging prescription, we construct age and 
metallicity histograms for the M31 globular cluster sam- 
ple, which are shown in Figure \7\ The upper left panel 
shows the age distribution of observed M31 globular clus- 
ters. A clear peak around 10 — 12 Gyr dominates the age 
structure. The majority of the globular cluster sample 
(~ 60%) has formal ages older than 9 Gyr. A surprising 
feature of the age structure in the observed M31 globular 
cluster system is a population of intermediate-age objects 
with formal ages between 6 and 8 Gyr. The mean un- 
certainty at those ages varies between 2 and 3.5 Gyr, so 
that it is unlikely that these objects have much older or 
younger ages. In particular, this feature is not reproduced 
by the overlap of old iso-age tracks, which is the result of 
increasing flux from hot horizontal-branch stars with in- 
creasing age, at metallicities below ~ —0.6 dex. Only six 
globular clusters are affected by the region where the age 

Note that the lBruzual fc Charloli i2003l) models do not in- 
clude predictions for stellar pop ulations with a given [a /Pe l 
ratio, contrary to the models of iThomas et all ll2003l 12004*). 
However, given the young ages, the sensitivity of the H/3 index 
to [a/Pe] variations is relatively weak, assuming the changes 
in chemical composition are correctly m odeled. Hence, t o first 
order the age-metallicity predictions of IBeaslev et alj ll2004fl 
should be relatively robus t. This is not the case for higher- 
order Balmer Lick indices ijPuzia et al.ll2004hl) . 



determination is ambiguous. For these clusters the code 
randomly assigns either old or young age values. 

After clipping globular cl usters with ages < 1 Gyr, a 
heteroscadastic KMM test* l)Ashman et al.lll994l) shows 
that a bimodal age distribution is more likely than a 
single-peak Gaussian distribution with a mean at ~ 11 
Gyr. We find a 0.1% likelihood for the single-peak dis- 
tribution and point out that a bimodal distribution is 
also more likely than a tri- modal fit. For the bimodal dis- 
tribution the KMM test gives peaks at ti = 10.9 ± 0.3 
and t2 — 6.5 ± 0.4 Gyr^ with dispersions 1.7 and 3.5 
Gyr, respectively. The integrity of the intermediate-age 
peak persists when we overplot a non-parametric proba- 
bility density estimate with its 90% confidence limits. We 
use a non-parani etric variable-width Epanechnikov kernel 
l(Silvermanlll986l) that is consistent with the mean age un- 
certainty of each individual globular cluster. 

In the left middle-row panel of Figure we present 
the metallicity distribution of studied M31 globular clus- 
ters. The distribution is clearly bimodal stretching from 
around —2.2 dex to values slightly above solar. A het- 
eroscedastic KMM-test gives peaks at — 1.66 ± 0.05 and 
-0.45 ± 0.04 dex with dispersions 0.23 and 0.29 dex, re- 
spectively. From an empirical calibration o f optical/near- 
IR colors with Milky Way globular clusters. lBarmbv et 
(2000) find peaks around -1.4 ± 0.05 and -0.6 ± 0.05 
dex^°. We note that our sample is biased towards metal- 
rich globular clusters which is a consequence of the sam- 
ple selection. The latter was driven by the observational 
setup and focused mainly on M31's d isk and bulge (see 
IPerrett et a'l]l2002t IBeaslev et al.ll2004 . Hence, the mean 
of each peak and its dispersion are subject to change for 
larger samples. In particular, the ratio between metal-poor 
and metal-rich globular clusters is not expected to be rep- 
resentative for the entire globular cluster system. This is 
at least partl y the reason for the offset between our spec- 
troscopic and lBarmbv et all 's photometric mean metallic- 
ity in the metal-poor peak. However, it is less likely to 
explain the ~ 3a offset for the metal-rich peak entirely, 
which might be a consequence of different absolute metal- 
licity scales of the used SSP models. 

In the upper right panel of Figure [3 we present the 
age-metallicity correlation plot. We find old globular clus- 
ters at all metallicities. For globular clusters with formal 
ages below ~ 8 Gyr, the evidence for metallicity bimodal- 
ity remains. A group of five clusters (B315-038, B321-046, 
B322-049, B327-053, and B380-313) with ages below - 5 
Gyr and metallicities around —1.8 dex appears as a very 

* Two Gaussian modes with independent mean values and 
dispersion are fit to the data, as opposed to a homoscedastic 
fit where two modes with independent mean values but same 
dispersions are adopted during the fit. 

^ The given errors are errors of the mean. Systematic errors 
are larger and of the order of 2-3 Gyr. Por a homoscedastic 
KMM test the peak values change only marginally within the 
uncertainties. The mean dispersion of both peaks is 2.4 Gyr. 

These metallicity values are on the Zinn-West scale, which 
is generally adopted throughout this work. 
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distinct feature in the plot. The other metal-rich young 
objects are clustered around [Z/H] « —0.6 to —0.2 dex. If 
the group of five young metal-poor globular clusters is ex- 
cluded, the overall age dispersion increases from amp « 1.3 
Gyr for objects with metallicities below —1.0 dex to crn„. ~ 
3.4 Gyr for metal-rich counterparts. The mean metallicity 
of intermediate-age globular clusters with ages between 3 
and 8 Gyr is —0.61 ± 0.11 dex with a dispersion of 0.48 
dex. Without two clusters from the young metal-poor pop- 
ulation (B321-046 and B380-313), the mean changes to 
— 0.49±0.08 dex and a dispersion 0.34 dex. There is tenta- 
tive evidence that this relatively metal-rich intermediate- 
age globular cluster sub-population might have a bimodal 
metallicity distribution in itself with peaks around —0.6 
and —0.2 dex. Such a multi-modality in the metallicity 
distribu tion of M31 globular clusters was already discov- 
ered bv iBarmbv et alJ l)2000(l in their photometric study 
(see their Fig. 19). It is also worthy of notice that the 
metal-poor peak of this intermediate-age sub-sample co- 
incides with the metal-rich peak of old globular clusters, 
while the metal-rich intermediate-age clusters appear to 
have ~ 0.4 dex higher mean metallicities than the aver- 
age metal-rich old globular cluster. However, these results 
need follow-up analyses with larger samples. 

Irrespective of the completeness and metallicity cov- 
erage of our sample, we do find subpopulations of old 
(~ 10— 12 Gyr), intermediate-age (~ 6 — 8 Gyr), and few 
young 2 Gyr) globular clusters in M31. Our results are 
fully consistent with the re cent spectro-photo metric study 
of 172 globular clusters bv l.Tiang et al.l l|2003l) . who find a 
significant population of intermediate-age and young glob- 
ular clusters (see their Fig. 6). Their results also indicate 
that metal-rich globular clusters are on average younger 
than metal-poor globular clusters. It would be presump- 
tions to speculate on the significance of this result ex parte, 
given our limited sample size (70 out of 460 ± 70 globular 
clusters) . With the current dataset we cannot conclusively 
estimate the fractions of these new types of intermediate- 
age and young globular clusters and assess whether these 
sub-populations are small components or fundamental in- 
gredients of the M31 globular cluster system. In any case, 
the distribution of globular cluster ages suggests that M31 
managed to form star clusters that survived as globular 
clusters, until a few Gyr ago. 

4.1.2. [a/Fe] Ratios 

The determination of [a/Fe] ratios yields very stable val- 
ues within a few iterations of our grid-interpolation rou- 
tine, which is the result of the small dependence of the 
[a/Fe] diagnostic grid on age and metallicity. The left 
bottom-row panel of Figure shows the distribution of 
[a/Fe] ratios, for which a single-peak distribution appears 
to be a good approximation. The mean of the single- 
peak distribution is 0.14 ± 0.04 dex^^ with a dispersion 
a = 0.32 dex. However, using the KMM test we find that 

Median: 0.18 dex with 0.15 dex semi-interquartile range. 
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Fig. 7. Ages, metallicities, and [a/Fe] ratios for M31 glob- 
ular clusters. Left panels show histograms for age, metal- 
licity, and [a/Fe] ratios for our sample globular clus- 
ters. Solid lines represent non-parametric probability den- 
sity estimates us ing a variable-width Epanechnikov kernel 
ijSilvermarJ 1 9861) with their 90% confidence limits (dashed 
lines). Right panels show correlations between the three 
parameters. Error bars are not shown for clarity; the mean 
error in age is Ai/i « 1/3, the mean uncertainty for metal- 
licity and [a/Fe] is ~ 0.2 dex. 



a two-Gaussian distribution is more likely (99.4%) than 
the single-peak distribution. The test yields a relatively 
narrow peak at [a/Fe] = 0.20 ± 0.02 dex with a dispersion 
(7 = 0.19 dex, and a relatively broad peak at —0.03 ±0.16 
with a factor 2.5 larger dispersion of cr = 0.50 dex. The 
mixture of these modes is six to one, where the super-solar 
[a/Fe] peak is more populated. 

We can now correlate [a/Fe] ratios of our sample glob- 
ular clusters with their ages and metallicities. The right 
middle-row and bottom-row panels of Figure \7\ illustrate 
the metallicity vs. [a/Fc] and age vs. [a/Fe] plots. We find 
no correlation in the metallicity vs. [a/Fe] plot. Metal- 
poor and metal-rich clusters have similar super-solar mean 
[a/Fe] ratios within the uncertainty of 0.1 — 0.2 dex. We 
find an increasing dispersion towards lower metallicities, 
which is mainly due to a decreasing [a/Fe] resolution at 
lower [Z/H]. The sample of the five clusters with ages 
~ 1 Gyr and metallicities around —1.6 dex shows super- 
solar [a/Fe] ratio +0.21 ±0.11. Metal-rich ([Z/H] > -1.0 
dex) globular clusters with ages below 8 Gyr show solar- 
type a/Fe ratios, i.e. [a/Fe] = 0.05 ± 0.05 dex. Their 
old counterparts, on the other hand, have a higher mean 
[a/Fe] — 0.23 ± 0.03. Hence, we find evidence for an age- 
[a/Fe] relation in the metal-rich globular cluster system. 
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An eye-catching feature of the age- [a/Fe] correlation 
plot (right bottom-row panel of Figured) is the decreasing 
scatter in [a/Fe] with decreasing age. Globular clusters 
with ages above 8 Gyr have a mean [a/Fe] = 0.18 ± 0.05 
and a dispersion of 0.37 dex^^. Younger objects have a 
smaller [a/Fe] ratio of 0.09 ± 0.09 and a dispersion of 0.21 
dex. The sub-sample of globular clusters with ages ^ 2 
Gyr is consistent with solar abundance ratios [a/Fe] — 
0.03 ± 0.08 dex and a cr = 0.18 dex. 



4.2. CN Enhancement 

In the following section, we investigate correlations be- 
tween ages, metallicities, and [a/Fe] ratios and other 
Lick indices for our sample of M31 globular clusters. 
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correlations of the CN index with globular cluster age, 
metallicity, and [a/Fe] ratio. 

Since CN2 has its blue bac kground passband d esigned 
to avoid the RS (see Fig. 3 in lPuzia et al]l2004b|) feature 
and has a calibration uncertainty of the transformation to 
the Lick system which is smaller than for CNi, we favor 
CN2 throughout the subsequent analysis and refer to it as 
the CN index. The CN index is a measure of the strength 
of the molecular CN band at 4150 A. 

In Figure IHl we plot CN vs. [MgFe]' along with SSP 
models for ages from 1 to 15 Gyr (isochrones are solid 
lines), metallicities —2.25 to -fO.67 dex (iso-metallicity 
tracks are dashed lines), and a [a/Fe] ratio of 0.3 dex. 
The model grid is completely degenerate in age, metallic- 
ity, and [a/Fe] and serves only as a reference line. In the 
plot we parameterize M31 globular clusters by their age: 
open circles mark old clusters with formal ages > 8 Gyr, 
open circles with a central plus sign are intermediate-age 
globular clusters with ages between 5 and 8 Gyr, and solid 
circles show clusters formally younger than 5 Gyr. We find 
a clear offset in the CN-[MgFe]' sequence between old and 
young M31 globular clusters. Young objects tend to have 
weak CN indices and fall below the model tracks, while 
old globular clusters show strong CN indices and populate 
the region above the model tracks. The CN index offset 
in the range [MgFe] '= - 1.8 A ([Z/H] between -2.3 and 
—0.6 dex) between these two groups is 0.13 ±0.02 A. The 
intermediate-age group falls in between the young and old 
sub-population. However, a ACN = 0.11 ± 0.03 with re- 
spect to young globular clusters makes it rather consistent 
with CN indices of old globular clusters. An interesting 
observation is the scarcity of globular clusters with small 
CN indices at [MgFe] > 2 A. The mechanism which drives 
the CN difference between young and old globular clus- 
ters may be a function of metallicity. Unfortunately, our 



sample is still too small to explore any metallicity depen- 
dence. 



4.2.1. Is the CN-ofFset a Consequence of Age, 
Metallicity, and/or [a/Fe] Differences? 

The CN index is mainly sensitive to the carbon and nitro- 
gen abundance. In relatively cool giants that dominate the 
integrated light of old stellar populations, most of the car- 
bon is locked in CO, so that the nitrogen abundance affects 
the CN index more directly. Main-s equence dwarfs do no t 
contribute to the CN absorption ijBurstein et al.l ir984j) . 
since molecules dissociate at higher log 5. Smaller CN ab- 
sorption is also consistent with higher effective tempera- 
tures, which are expected in young a nd /or metal-po o r stel- 
lar populations. The SSP models of iThomas et alJ l)2003l 
I2OO4) predict an offset of ACN « 0.02 - 0.04 mag be- 
tween 13 Gyr and 3 Gyr^^ for stellar populations with 
metallicities between [Z/H] « —1.35 and 0.0, and solar- 
type abundance ratios. Since the models were calculated 
for defined [a/Fe] ratios, we can test the influence of a- 
element enhancement on the CN index. The maximum 
CN index offset between stellar populations with [a/Fe] 
0.5 and 0.0 dex, for ages up to 13 Gyr and metallicities up 
to solar, is ~ 0.02 mag. For super-solar metallicities the 
difference does not exceed 0.025 mag. Age and metallicity 
are therefore excluded as primary parameters that drive 
the CN offset. 

Since we find evidence for an age- [a/Fe] relation, we 
have reason to assume that the CN offset is produced 
by varying [a/Fe] ratios. For an [a/Fe] ratio of 0.2 dex 
the models predict very similar maximum CN offsets of 
~ 0.045 mag between 3 and 13 Gyr for solar and super- 
solar metallicites, still a factor 3 too small for the mea- 
sured offset. The models predict smaller offsets for smaller 
metallicities. Going to highly super-solar [a/Fe] ratios 
^ 0.5 dex, we find CN offsets of the order of 0.1 mag. 
However, such high a-element enhancements are incon- 
sistent with the measured [a/Fe] ratios of the young sub- 
sample. In summary, any combination of differences in age, 
metallicity, and [a/Fe] ratio between young and old M31 
globular clusters cannot explain the offset in their mean 
CN index. 



4.2.2. Is the CN-offset a Consequence of Carbon 
and/or Nitrogen Abundance Differences? 

Because iThomas et all l(2003j) provide index predictions 
for stellar populations with an enhanced nitrogen abun- 
dance, we can compute a zeroth-order estimate of the in- 
crease in nitrogen abundance that is consistent with the 
observed CN index offset. Table 13 summarizes offsets in 
the CN index as a response to a factor three increased 
nitrogen abundance for different ages, metallicities, and 
[a/Fe] ratios. Note that the abundances of other species 



These values exclude two outliers at [a/Fe] 
dex. 



4.0 and 1.0 



The mean age of the sub-population of globular clusters 
with ages below 5 Gyr is 2.5 ± 1 Gyr. 
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Fig. 8. CN2, G4300, C24668, and Ca4227 vs. [MgFe]' for M31 (drcZes), Milky Way {squares), M81 (diamonds), 
LMC (stars), and Sculptor-group spiral galaxy globular clusters (pentagons). The individual datasets are described 
in Section ^21 The M31 symbols are parameterized by the globular cluster age: filled circles show globular clusters 
younger than 5 Gyr, circles with a central plus sign are intermediate-age clusters with ages between 5 and 8 Gyr, while 
open circles show globular clusters older than 8 Gyr. Note the offsets in CN and Ca4227 indices between old (open 
circles) and young (filled circles) M31 globular clusters which is most likely a result of a factor three relative nitrogen 
enhancement of the old cluster system. 



remain unchanged. While other indices do not respond sig- 
nificantly, the increase in nitrogen abundance is expecte d 
to decrease the Ca4227 index (see iTrinicco fc Beiilll99,ql) . 
Indeed, we observe such a decrease which is illustrated 
in Figure |H| However, we refrain from a more quantita- 
tive analysis due to the poor interna l calibration quality 
of this index (see IPuzia et al.ll2002bl) . For the mean age 
~ 13 Gyr, a metallicity ~ —0.4 dex, and [a/Fe] ratio ^ 0.2 
dex, we find a predicted CN offset of the order ~ 0.07 mag, 
consistent with a factor three increase in nitrogen abun- 



dance. Together with the CN offset resulting from the age 
difference and under the assumption that the young glob- 
ular cluster sub-population is not enhanced in nitrogen, 
the predicted CN offset amounts to ^ 0.09 — 0.11 mag, 
which is in the right ballpark as the measured value. In 
fact, the nitrogen enhancement appears to be even larger 
than a factor of three. 

Since the Lick system provides indices which are pri- 
marily sensitive to carbon abundance, we attempt to ad- 
dress the question of whether a variation in carbon abun- 
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Table 3. Index response of the CN index to a factor three 
increase in nitrogen abundance. The valu es were calcu- 
lated with model predictions taken from omas et al.l 

(Eool. 
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0.030 
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dance can explain the CN-index offset between young and 
old M31 globular clusters. For this purpose we plot G4300 
and C24668 vs. [MgFe]' in the middle and right panel 
of Figure |H1 G4300 measures the strength of the G-band 
which is dominated by CH-band absorption and a few Fe 
and Cr features; C24668, on the other hand, measures the 
C2 Swan band at 4700 A, but its index passbands also in- 
clude many weak Fe, Ti, Cr, and Sc lines. Both indi ces are 
prim arily sensitive to carbon abundance iTripicco fc Belli 
To first order we do not see any significant off- 
set between young and old globular clusters in the CN 
and C24668 index plot, although there are indications for 
differences in G-band strength at very low [MgFe]' (i.e. 
metallicities) . We therefore conclude that a difference in 
carbon abundance is not the driving factor for the ob- 
served CN offset. 

A more detailed look at the red background passband 
of the CN index reveals that the index is contaminated 
by two iron features at AA4260 and 4272 A while the blue 
background passband has a weak Mn i feature located at 
its blue edge at A4084 A. Although this Mn feature is fairly 
weak, changes in Mn abundance might have a significant 
influence on the CN index. The blue background passband 
is relatively narrow (12.5 A) and located ~ 50 A away from 
the main CN feature, and small changes in Mn abundance 
can therefore translate in a large leverage of the pseudo- 
continuum flux inside the feature passband. Hence, com- 
plex, non-solar abundance ratio variations might influence 
the CN index strength, as well. However, a quantitative 
analysis of these effects is beyond the scope of this work 
as detailed modeling is not available at the time. 

We conclude that the offset in CN index between young 
and old M31 globular clusters is likely due to a nitrogen 
enhancement of old globular clusters by a factor three or 
more, compared to the young globular clust er system. Our 
result s are consistent with the s tudies of iLi fc BursteinI 
l|2003(l and iBurstein etaP l|2004|) who find a systemati- 
cally higher nitrogen abundance in M31 globular clusters 
as compared to Galactic counterparts^^. 



4.3. Comparison with Globular Clusters in other Spiral 
and Dwarf Galaxies 

In this Section we present a comparison of globular clus- 
ters in M31, Milky Way, M81 (Sa), LMC, and Sculptor- 
group spiral (Sa-b) galaxies. To avoid uncontrolled biases 
due to varying instrumental setups and spectral charac- 
teristics (i.e. resolution, sky-subtraction, flux-calibration, 
etc.), we strictly compare parameters derived from in- 
dex measurements that were obtained in the same in- 
dex system. For conformity with SSP mo del predictions 
we ch oose index passbands defined by IWorthev et alJ 
1I1994I) and for higher- order Balmer lines the passbands of 
IWorthev fc Ottavianil (|l997l) . The size of systematic un- 
certainties between different datasets is set by the cali- 
bration quality of the transformation to the Lick system. 
For this reason we collect only high-quality data from 
the literature with very good index calibrations to keep 
the systen iatics small. For the Milky Way sample we use 
data from IPuzia et alJ l|2002bl) i5_ 

and colle ct index mea- 
surements for LMC globular clusters from i Beaslev et all 
i2002). M81 globular cluster Lick index measurements 
were taken from ISchroder et afJ lj2002l) . For a handful 
of globular clusters in the Sculptor-group spiral galaxies 
NGC 55, 247, 253, and 300, Olscn ct al. (200^ measured 
Lick indices. A list of index measurements was kindly pro- 
vided to us by Knut Olsen. Due to small sample sizes we 
merge the M81 and Sculptor-group globular clusters into 
one sample with ten globular clusters in total. We warn 
however that the individual globular cluster samples can 
have very different age, metallicity, and abundance ratio 
distributions. Diagnostic diagrams for all globular clusters 
are shown in Figures |S1 and |5| 

4.3.1. Ages and Metallicities 

Figure El shows M31, Milky Way, M81, Sculptor-group spi- 
ral galaxy, and LMC globular clusters in three different 
age-metallicity diagnostic plots, as circles, squares, dia- 
monds, pentagons, and stars, respectively. After extract- 
ing ages, metallicities, and [a/Fe] ratios for globular clus- 
ters in other galaxies in exactly the same way as described 
earlier for M31, we find several interesting systematic dif- 
ferences between the observed globular cluster systems. 
Age, metallicity, and [a/Fe] distributions are shown in 
Figure UOl Both the M31 and Milky Way globular clus- 
ter systems are dominated by objects reaching oldest ages 
at ~ 11 — 13 Gyr. In addition to old clusters, the M31 sys- 
tem hosts populations of intermediate-age (6-8 Gyr) and 
young {<, 2 Gyr) globular clusters, which are not observed 
in the Milky Way sample. 



^'^ We note that the metallicity coverage of their globular clus- 
ter sample is not consistent between M31 and Galactic globular 
clusters. This point is discussed further below in the Discussion 
section. 



Although we are aware of the^Cohe n et al.l 1^9^ data, we 
do not use these in the following analysis because of the lack of 
H(5a index measurements, which might introduce systematics 
in our analysis. However, the available index measurements 
are included in following diagnostic plots and marked as dotted 
squares. In general, they agree very well with the m ean locus of 
the other globular cluster data taken from lPuzia et ah (.2002bl) . 
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Fig. 9. Comparison of age-metallicity diagnostic diagrams panels a-c) for M31 {solid circles), M81 (diamonds) , Milky 
Way (solid squares), and LMC globular clusters (stors), as well as globular clusters in Sculptor-group spiral galaxies 
(pentagons). The panels show diagnostic plots using different B ai mer indices. H74, H/3 , and US a- Also shown is the 
a/Fe diagnostic grid panel d). Milky Way globular clusters are from lPuzia et alJ l)2002tJ) . while LMC globular clusters 
were taken fr qm.,Bcaslev ct al. (200^). Squares with a central dot indicat e measurements of Galactic globular clusters 
performed by 
performed by 



Cohen et alJ iflQQSl) a nd cal ibrated bvlBeaslev et all ((20041) . Lick indices for M81 globular clusters were 
Schroder et alJ l)2002l) . and lOlsen et alJ ^|2004^ measured indices for globular clusters in Sculptor-group 
spiral galaxies. Note that USa measurements are not available in the Cohen et al. dataset. Most of the M81 data 
have no higher-order Balmer index measurements. Also plotted are index measurements for the nuclear diffuse light of 
the Milky Way and M31. A large open s quare marks the index measurements for the Galactic Bulge light in Baade's 
Window, taken from P uzia et alJ l|2002lJl . and two crossed open circles indicate measurements for the diffuse light of 
the M31 bulge at two different galactocentric radii (labeled as zO and zl, see Table^). Note, that all data and models 
use the same index passband system. SSP model predictions are as in Figure El 



The mean age of the Milky Way sample is 12.0 ± 1.3 
Gyr^^ which is in good agreement with previous CMD re- 



This uncertainty is the formal statistical uncertainty de- 
rived from our approximation routine, while the systematic 



suits (11.5 ± 1.3 Gvr. lChabover et al.lll998|) . The metallic- 



uncertainty derived from the comparison of results from dif- 
ferent diagnostic plots is ~ 3.5 Gyr. The median value is 12.1 
Gyr. 
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ities covered by M31 and Milky Way globular clusters are 
similar, which again underlines the importance of careful 
sample selection when different globular cluster systems 
are compared (see Sect. \2A\i . The metallicity distribution 
of the Milky Way sample is bimodal with peaks around 
— 1.2 and —0.3 dex, although these estimates are uncertain 
due to the small sample size of each peak. 

LMC globular clusters span a wide range in age from 
objects as young as '--^ 1 Gyr to old clusters resembling 
ages of old Galactic counte rparts. This is consi stent with 
previous findings (see also iBeaslev et al]l2002l and ref- 
erences therein) . The metallicity distribution of the LMC 
sample is consistent with a broad single peak with a mean 
[Z/H] = -0.95 ±0.09 dex^"^ and a dispersion a = 0.42 dex. 

We combine the samples of M81 and Sculptor-group 
globular clusters because of the small number of stud- 
ied objects in each galaxy. The combined age distribution 
shows peaks around 12 and 8 Gyr with a clear indication 
for intermediate-age globular clusters, which is similar to 
the M31 rather than the Milky Way globular cluster age 
distribution. However, very young objects are missing in 
this sample. The metallicity distribution is peaked around 
—0.6 dex with only one cluster at roughly solar metallic- 
ity and one metal-poor candidate at —1.5 dex. We warn 
that the distributions are subject to change once larger 
and more representative samples are available. 

Although the samples of observed globular clusters in 
all galaxies are still small, all cover the brightest fraction 
of each globular cluster system. Keeping this in mind, we 
can draw the following conclusions. The age structure of 
each globular cluster system suggests that luminous glob- 
ular clusters experienced significantly different formation 
histories in all three galaxies. While the question about 
the presence of intermediate-age globular clusters in the 
Milky Way will certainly be answered by our ongoing spec- 
troscopic survey of Local Group globular clusters, photo- 
metric CMD studies suggest that a small fraction of metal- 
rich Galactic globular clust ers might be younger than ~ 9 
Gyr l|Salaris fc Weis jl2002l) . However, among the studied 
bright clusters we find no Galactic counterparts for the 
intermediate-age and/or young globular clusters found in 
M31, M81, LMC, and Sculptor-group spiral galaxies. From 
Figure [TUl it seems that the age structure of the M31 clus- 
ter system is a composite of globular clusters in all other 
studied galaxies. This aspect might suggest that M31 glob- 
ular clusters with ages ^ 9 Gyr were accreted from satel- 
lite galaxies. 

4.3.2. [a/Fe] Ratios 

The average [a/Fe] ratios of LMC and Milky Way glob- 
ular clusters are ~ 0.1 — 0.2 dex higher than the mean 
a-element enhancement of M31 globular clusters (see 
also Fig. inii)- The mean a-element enhancement of the 
Galactic sample is clearly super-solar with a mean 0.35 ± 




15 10 5 -2 -1 -1-0.5 0.5 1 
t [Gyr] [Z/H] [a/Fe] 

Fig. 10. Age, metallicity, and [a/Fe] distributions for 
globular clusters in M31, Milky Way, LMC, M81, and 
Sculptor-group spiral galaxies. Due to small sample sizes, 
the M81 and Sculptor-group samples were merged. Solid 
lines are non-parametric probability density estimates 
with their 90% bootstraped con fidence limits shown as 
dashed fines (see ISilvermanlllQS^ for details). 



0.04 dex^® and a dispersion cr = 0.13 dex. The mean [a/Fe] 
ratio of the LMC cluster sample is 0.32±0.04 dex and has a 
dispersion of 0.17 dex. The combined M81/Sculptor sam- 
ple has a mean [a/Fe] = —0.05 ± 0.08 and a fairly broad 
distribution with a dispersion of a = 0.25 dex. 

Yet, all three globular cluster systems show super- 
solar [a/Fe] ratios, which imp lies enrichment d ominated 
by type-II supernovae (e.g. iMatteiiccil Il994l) . Super- 
solar mean a-element enhancements around ~ 0.2 — 
0.35 dex are observed in high-resolution spectra of sin- 
gle stars in Galacti c halo and bulge g l obular clusters 
barbuv et all Il999l: ICohen et all Il999t ICarretta et alJ 
l200lLl2004HOrigfia et alJl2002l:IOrigfia fc Richll2004|) . This 
is in good agreement with the values derived here and 
can be conside red as a zero-point cal i bration of our [a/Fe] 
scale (see also iMaraston et al.ll2003t iThomas et al.ll2003l 
\2()04\ . However, there are Galactic globular clusters, such 
as Pal 12, which show peculi ar abundance signatures with 
roughly solar [a/Fe] ratios ijCoheiJ l2004l) . Judging from 



'^'^ The median of the distribution is —0.89 dex. 



The median value is 0.37 dex. 
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their kinematics, these objects appear to be stripped from 
nearby satellite galaxies. 

Compared to the mean [a/Fe] of LMC clusters, the 
mean [a/Fe] ratio of intermediate-age M31 globular clus- 
ters is ~ 0.1 — 0.2 dex lower, and suggests that accre- 
tion of LMC-type globular cluster systems to build up 
the intermediate-age globular cluster population in M31 
is unlikely. This low [a/Fe] is rather consistent with the 
mean value found for the M81/Sculptor cluster sample. It 
is also consistent with the abundance patter n in thin/thick 
disk stars in the Milky Way (e.g. lEdvardsson et alJll993t 
lFuhrmannlll998l) . 



4.3.3. Carbon and Nitrogen Abundance 

Beginning with the study of iBurstein et al.l l)l984|) . glob- 
ular clusters in M31 were found to be unique in individ- 
ual line indices compared to their Galactic counterparts. 
Enhanced CN-band absorption in M31 clus ters was later 
found to be due to nitrogen enhancemen t dPonde r et alJ 
Il99at iLi fc BursteinlEm iBurstein et al.l Eo04). How do 
the CN indices of M31 globular clusters compare to other 
globular cluster systems? In Figure [S] we plot all other 
globular clusters together with M31 globular clusters, pa- 
rameterized by their age. 

Old M31 globular clusters appear on average CN- 
enhanced compared to the other globular cluster systems. 
Intermediate-age clusters, on the other hand, compare well 
in their mean CN index with the other cluster systems. We 
find a striking coincidence in CN strength between LMC 
and young M31 globular clusters. The youngest M31 clus- 
ters populate the same locus in the CN vs. [MgFe]' plot 
as most LMC globular clusters. Compared to SSP models, 
which predict CN strength for solar-type abundance pat- 
tern, this suggests that both LMC and young M31 globu- 
lar clusters are underabundant in carbon and/or nitrogen. 
Since no systematic difference in the C24668 and G4300 
indices, which are most sensitive to carbon abundance, 
is found between LMC and young M31 globular clusters, 
the CN offset suggest a difference in nitrogen abundance. 
Consistent with this picture is the systematic Ca4227 off- 
set between old M31 globular clusters, on the one hand, 
and LMC and young M31 globular clusters, on the other 
hand (see Sect. 14.2(1 . 

4.4. Kinematics 

In a recent radial-velocity study o f the M31 globular clus- 
ter system, iMorrison et al.l l)2004() found a significant sub- 
population of globular clusters showing thin-disk kine- 
matics with metallicities down to —2.0 dex. The pres- 
ence of such a distinct disk globular cluster sub-system 
with low-metallicity members implies an early formation 
epoch of M31's thin disk. We correlate the kinematics of 
IMorrison et alJ with our ages, metallicities, and [a/Fe] ra- 
tios. A cut in normalized residual radial veloci tv at 1 0.75 1 
km s^^ relative to the thin-disk model of iMorrison et all 
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F ig. 11. Residual radial velocities relative to a disk model 
of IMorrison et all l)2004j) as a function of age, metallicity, 
and [a/Fe] ratio for M31 globular clusters. Error bars are 
not shown for clarity; the mean error in age is At/t « 1/3, 
the mean uncertainty in metallicity and [a/Fe] is ^ 0.2 
dex. The typical error in radial velocity is < 12 km s~^. 



I|2004|) selects globular clusters which are likely members 
of the thin-disk population, as suggested by the authors. 
This selection leaves 18 objects with thin-disk kinematics 
in the sample. 

This sub-sample covers ages between one and 11.9 Gyr. 
It has a mean age 8.1 ± 0.9 Gyr and a dispersion a = 
3.6 Gyr. The average metallicity is [Z/H] = -0.78 ± 0.16 
dex with a = 0.68 dex; its mean [a/Fe] is slightly super- 
solar at -1-0.12 ± 0.07 dex with a sample dispersion of 0.30 
dex. We confirm that thin-disk globular clusters span a 
wide range in metallici tv from about — 1.8 dex up to solar 
values, as reported bv iMorrison et al.1 (|2004) . 

If globular clusters with residuals smaller than |0.25| 
km s~^ are selected, six objects remain in the sample 
and all but one^^ h ave in termediate to young ages ^ 8 
Gyr. iBeaslev et al.l ()2004|) found also three young globu- 
lar clusters with clear thin-disk kinematics. Their median 
metallicity decreases to —1.02 ± 0.31 dex with a higher 
metallicity dispersion of 0.75 dex and their median [a/Fe] 
ratio increases to 0.26 ± 0.05 dex (u = 0.12 dex). This 
suggests that M31 globular clusters with clear thin-disk 
kinematics are members of the intermediate-age cluster 
sub-population. Although their metallicities are inconsis- 
tent with typical abundances found in the Milky Way's 
thin-disk stellar population, their [a/Fe] ratios are in ac- 
cord with those of Galactic bulge globular clusters (see 



This cluster (B171) has a formal age of 10.8 Gyr and might 
be a halo interloper. 
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Sect. 14. 3^ . It would be interesting to obtain a larger sam- 
ple of thin-disk globular clusters in M31 (1) to test the 
maximum age of this globular cluster sub-population, and 
(2) to check whether or not their chemical abundances are 
in contrast with disk/bulge globular clusters in the Milky 
Way. 

In Figure ^2 we plot the residual radial velocity, 
which was calc ulated relative to the thin disk model of 
iMorrison et al.l as a function of globular cluster age, metal- 
licity, and [a/Fe] ratio. There is a clear bimodality in the 
residuum distribution and we split the sample into clus- 
ters that fall in the hatched region which indicates glob- 
ular clusters with normalized disk-model residuals < |20| 
km s~^, and globular clusters with larger residuals indi- 
cating halo membership. We consider the sub-population 
with small residuals as a thin/thick disk component, being 
aware of potential contamination by halo objects. About 
70% of our globular cluster sample falls into the thin/thick 
disk category, excluding only halo members. Based on the 
number counts of the halo and thin/thick-disk sample and 
assuming a Gaussian distribution of velocities around the 
mean of the thin-disk model with dispersions 30 and 150 
km ss^^, for the disk and halo population, respectively, we 
expect 2±2 halo interlopers in our thin/thick-disk sample 
and less than one genuine disk globular cluster in the halo 
sample. 

The most striking feature in the upper panel of 
Figure^]is that most globular clusters of the halo popula- 
tion have ages older than 9 Gyr, with a median age of 10.4 
Gyr. Only two halo globular clusters have intermediate 
ages and we find one extremely young cluster (B324-S51) 
that is inconsistent with disk kinematics. The thin/thick- 
disk population, on the other hand, has a median age of 
9.4 Gyr. The halo sub-sample shows a clearly bimodal dis- 
tribution of metallicities with peaks around —1.6 and —0.5 
dex. 

It is interest ing that the velocity dispersion around the 
Morrison et al.l model in the thin/thick-disk sample de- 
creases towards younger globular clusters; from |0.47| km 
s~^ for objects with ages > 8 Gyr to 0.38 km s^^ for 
globular clusters younger than 8 Gyr. 

The distribution of [a/Fe] ratios for halo and 
thin/thick-disk globular clusters is shown in the lower 
panel of FieurelTTI We find an average [a/Fe] — 0.17±0.07 
dex with a dispersion of 0.33 dex, for the halo sample. 
The thin/thick-disk sample has a similar mean [a/Fe] at 
0.14 ± 0.01 dex and a dispersion of 0.32 dex. 

We find a mild age- [a/Fe] correlation in the thin/thick- 
disk sample of the form [a/Fe] = (0.013 ± 0.011) • t + 
(0.076±0.122), with a rms of 0.12 dex for a weighted least- 
square fit, where t is given in Gyr. We compute weights for 
individual globular clusters from an average of upper and 
lower 1(7 uncertainties. Normally, these errors are asym- 
metric due to the skewness of diagnostic grids. No signif- 
icant correlation between [a/Fe] and metallicity is found 
for thin/thick-disk globular clusters. 



4.5. Structural Parameters 

In the following section we study the correlations between 
globular clus ter core ra di, Rc, half-light radii, Rhi, tidal 
radu, Rt fsee lKinsll962ll . and the projected galactocentric 
distance, rproj with previously derived ages, metallicities, 
and [a/Fe] ratios. We collect data on structural param- 
eters for M31 globular clusters from iBarmbv fc Huchral 
l|200l|) . for Milky Way globular clusters from the 2003 
update of the McMas tcr catalog^" ( IIa rrI^ll996l) . and for 
LMC clusters from lMackev fc Gilmorel l)2003l) . 

Since for extragalactic globular clusters only projected 
galactocentric distances are available, we need to project 
the 3-D positions of the Galactic globular cluster system 
on a 2-D surface on the sky for a fair comparison between 
Galactic and extragalactic cluster systems. Fortunately, 
the McMaster catalog provides the 3-D positions of indi- 
vidual Galactic globular clusters in the reference system 
of the Sun, so that the conversion to a projected galac- 
tocentric distance is easy. We simulate the view towards 
the Milky Way globular cluster system similar to that to- 
wards M31, by applying a similar pitch angle of 77° for 
the Milky Way disk and the entire globular cluster system 
and randomly choosing a line of sight tilted with respect 
to the Earth-Galactic center axis. 

In Figure 1121 we compare structural parameters with 
globular cluster ages, metallicities, and [a/Fe] ratios. We 
also plot the structural parameters as a function of pro- 
jected galactocentric radius in the very right column of 
Figure [T^ For each sub-panel we compute the Spearman 
rank correlation coefficient for the M31 globular cluster 
sample. The number is given in a corner of each panel (1 
indicates perfect correlation, —1 anti-correlation). 

In general, the parameter space coverage between the 
three globular cluster systems is in good agreement. We 
find correlations in the Rhi and Rt vs. [Z/H] panels, where 
the Spearman correlation coefficient is around —0.7 with a 
statistical significance for a correlation of > 99.7% in both 
cases. Since both the half-light radius Rm and the tidal 
radius Rt appear to be weakly correlated with galacto- 
centric radius Tpioj , at least part of their correlation with 
metallicity might be due to a correlation of metallicity 
with galactocentric radius (reflecting the galactocentric 
distance dichotomy between met al-rich disk and metal- 
poor halo globula r clusters; sec 'Cher noff fc Diorgovskil 
Il989t Ivan den Bergh,1991,: .Cote. 19991 This is consistent 
with the correlation in the rproj vs. [Z/H] plot in FigurelT^ 
However, the sample selection might produce spurious cor- 
relations and it is important to check these relations with 
larger samples, in particular because the Galactic sample 
does not show such a strong correlation of R^i and Rt 
with metallicity. 

We point out that similar size differences between 
metal-poor and n ietal-rich globular clusters were found 
in other galaxies llKundu fc Whitmorelll998t IPuzia et alJ 
hflflflt iT.arsen et a,1.ll200lh ~ possible explanation for this 

http:/ /physun. physics. mcmaster.ca/~harris/mwgc.dat 
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Fig. 12. Structural parameters, Rc (core radius), i?hi (half-light radius), and Rt (tidal radius), as a function of age, 
metallicity, and [a/Fe] ratio for globul ar clusters in M31 (circles). Milky Way (squares), and LMC (stars). Data on 
structural parameters \ vere taken fro m |Barmbv & Huchral (j^QOl) for M31 globu lar clusters, from th e 2003 update of 
the McMaster catalog l|Harrislll996 ^ for Milky Way globular clusters, and from iMackev &: GilmoT^ l)2003^ for LMC 
clusters. Information on age and m etallicity for Galactic globular clusters is parameterized by their photo metric (open 
square s, data taken from Salaris fc Weissll2002|) and spectroscopic origin (solid squares, data taken from IPuzia et alJ 
l2002bt) . The structural parameters are shown as a function of projected galactocentric radius, rproj, in the right-column 
panels. The bottom row of panels shows the projected galactocentric radius as a function of age, metallicity, and [a/Fe] 
ratio. In a corner of each panel the Spearman rank correlation coefficient is given. 



difference might be the combined effects of mass segre- 
gation and the dependence of main-sequence lifetimes on 
metallicity, under the assumption of similar distribution 
of central potentials and half- mass radii (Jordan 2004). 
However, the size difference might be also entirely driven 
by projection effects which are due to a systematic differ- 
ence in dy namical evolution of glo bular clusters on differ- 
ent orbits l|Larsen fc BrodielEooi . If true, size differences 
should disappear at large galactocentric radii, where dy- 
namical effects similarly affect metal-poor and metal-rich 
globular clusters. Intriguingly, Perrett (2001, PhD thesis) 
found an increasing size of metal-rich globular clusters 
with galactocentric radius, approaching the mean size of 
metal-poor globular clusters at large galactocentric dis- 



tances. No galactocentric distance-size relation was found 
in the metal-poor globular cluster sub-population. Size in- 
formation on more metal-rich globular clusters at large 
galactocentric radii is needed to establish whether the 
size difference between metal-poor and metal-rich globular 
clusters is due to a dynamical effect or reflects differences 
in globular cluster formation. 

It is well known that the globular cluster popula- 
tions of M31, LMC, and the Milky Way faU on the same 
" fundamental plane" , which for globular clusters is de- 
fined by velocity d ispersion, surface brightness, and radius 
l)Diorgovskilll995|) . Together with the similar correlations 
between the structural parameters of globular clusters in 
M31, Milky Way and LMC with their ages, metallicites, 
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and abundance ratios, this suggests that the formation 
and evolution must have been similar for the studied clus- 
ters. 

4.6. Comparison with the Integrated Bulge Light of 
M31 and the Milky Way 

Early speculations bv lBaadj l)l944) claimed that the stel- 
lar populations in Galactic globular clusters would resem- 
ble those in the Milky Way bulge and the nucleus of M31. 
This simplistic picture was revised soon after, when a 
more complete census of stellar populations in Galactic 
globular clusters was available. Only the most metal-rich 
globular clusters appear to host stellar populations simi- 
lar to those in spiral bulges and early-type galaxies (e.g. 
iMaraston et alJl2003(l . 

In the following section we compare the integrated- 
light spectra of the M31 nucleus, obtained i n this study 
(see Se ct.O, the Milky Way bulge (taken from lPuzia et al.l 
and metal-rich globular cluster populations of 
both galaxies. Figure shows the location of the two M31 
bulge spectra (as crossed open circles) in age/metallicity 
and [a/Fe] diagnostic grids for the two different galac- 
tocentric radii (see Fig. and Tab. QJ. Using the same 
technique as for the globular cluster spectra we derive 
luminosity-weighted ages, metallicities, and [a/Fe] ratios. 
The values are summarized in Table0] Since the inner M31 
field and the Milky Way field were taken at similar mean 
galactocentric radii ('^ 0.6 kpc), their stellar populations 
can be directly compared, assuming that no significant dif- 
ferences in population gradients are present between the 
two galaxies. 

Given the statistical and systematic errors, the derived 
mean ages, metallicities, and [a/Fe] ratios suggest that the 
stellar populations in both spiral bulges are similar. Both 
bulges harbor old stellar populations, with metallicities 
close to solar values. There is a tendency for the M31 
bulge to have a slightly higher metallicity. Both stellar 
populations have clearly super-solar a/Fe ratios, although 
the outer M31 field seems to have a higher [a/Fe] value. 

There is no systematic difference in CN index between 
metal-rich old M31 globular clusters, which are enhanced 
in nitrogen (see Sect. 14. 2|) . and the corresponding M31 
bulge fields (see Fig.lSj). However, compared to metal-rich 
Galactic globular clusters and the Galactic bulge, the CN 
index of the M31 bulge is offset by ^ 0.04 and ^ 0.1 
mag, respectively, and suggest an enhacement in carbon 
and/or nitrogen abundance. The C24668 index for the 
central M31 field (zO) is significantly higher with respect 
to the mean M31 and Galactic globular clusters at simi- 
lar [MgFe]'. This is not the case for the outer field (zl). 
However, this C24668 offset is less pronounced when the 
central M31 field is compared to the Milky Way bulge. In 
combination with the similar G4300 indices for metal-rich 
globular clusters and the bulges of M31 (both sectors) 
and the Milky Way, this suggests that the central M31 
field (zO) is enhanced in carbon and nitrogen, compared to 



Galactic bulge and metal-rich Galactic globular clusters. 
The outer M31 field is likely to be enhanced in nitrogen 
only. The reduced Ca4227 indices of both M31 fields with 
respect to metal-rich Milky Way globular clusters and 
Galactic bulge underline these results. Furthermore, we 
find an increase in CN index with increasing galactocen- 
tric radius for the M31 bulge. This suggests an increasing 
nitrogen enhancement, which in the case of the field stel- 
lar population in M31 increases with radius. Such differ- 
ences were already reported bv iMorgan fc Mavalll l|l95 /t) 
who discovered that the integrated spectrum of M31's 
nucleus is dominated by cyanogen absorption features. 
Re cently, radial NH grad ients in the M31 bulge were found 
bv lDavidge et alJ i|199(TI) . 

In summary, the bulge of M31 appears to be metal- 
rich, a-element enhanced and old, similar to the Milky 
Way bulge. However, compared to the Galactic stellar 
populations there are significant differences in chemical 
composition. We find indications for an enhancement in 
carbon and nitrogen in the central region of M31, while 
the integrated light of the outer stellar populations is 
consistent with a nitrogen enhancement only. It is well 
known that stellar populations in early-type galaxies show 
super-solar metalliciti es and a-element abundance ratios 
l|Worthev et alJ Il992^ . In this respect, our comparison 
shows that metal-rich M31 globular clusters and early- 
type galaxies are more alike than stellar populations in 
early-type galaxies and metal-rich Galactic globular clus- 
ters. 

5. Discussion 

5.1. Stellar Streams and Merging 

The existence of the recently disc overed thin-disk pop- 
ulation of M31 globular clusters (" Morrison et all l2004l) 
is a strong argument for an undisturbed disk evolution. 
Depending on the exact kinematical definition of the thin- 
disk population, we find that most, but not all, mem- 
bers of the thin-disk globular cluster population have ages 
younger than ^ 10 Gyr and metallicities —2.0 ^ [Z/H] ^ 
0.0 (see Sect. 14. 4| . The presence of old thin-disk globular 
clusters suggests an early formation and smooth kinemat- 
ical evolution of the thin-disk globular-cluster population. 
In particular, it implies in-situ formed star clusters. It 
would be therefore of great interest to determine the max- 
imum age in a larger dataset of thin-disk globular clusters 
and to test how early M31's thin disk formed. 

Recent very deep HST observations, pointed towards 
the halo of M31, made the age and metallicity distribu- 
tions of the targeted stellar populations directly accessible. 
These CMD studies revealed a dominant metal-rich stellar 
popul ation with a surprisingly high metallicity, fZ/H] » 
-0.6 ^Holland et al."l996' 'Saraiedin ijg Van Duvnjl200ll: 
iBellazzini ct al. 2003: Rich et al. 20M see also Fig. [B 
and intermediate ages ~ 6-8 Gvr JBrown et al J 1200.1) . 
There are indications for an even more metal-rich stel- 
lar population [Z/H] ^ —0.2) which is irregularly dis- 
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Table 4. Ages, metallicities, and [a/Fe] ratios of the integrated light for three fields in the bulge of M31 and the Milky 
Way. The spectra of all fields were taken at different median galactocentric distances, rmcdian, given in the second 
column. All measurements were performed in the same way as for the globular cluster data and are given with their 
statistical la uncertainties. 



galaxy 


'"median [kpc] 


age [Gyr] 


[Z/H] 


[a/Fe] 


MW 


0.55 


9.4 ± 1.4 


-0.22 ±0.06 


0.35 ± 0.07 


M31 (zO) 


0.60 


14.5 ± 1.3 


-0.06 ±0.04 


0.38 ± 0.07 


M31 (zl) 


2.48 


14.9 ± 2.2 


-0.22 ±0.12 


0.65 ±0.22 



tribti t ed throughout M31's h alo l|Saraiedini fc Van Duvnj 
l200lt iBellazzini et al.ll2003(l . The case of a metal-rich 
intermediate-age stellar population in the halo of a spi- 
ral galaxy is generally attributed to an accretion or merg- 
ing event that led to the incorporation of external stellar 
populations and/or to infall of fresh material which fu- 
eled new star formation. Depending on the exact merging 
kinematics, tidally disrupted satellites remain in the halo 
for several Gyr, until their central regions eventually spi- 
ral inwards as a result of dynamical friction. Support for 
this picture comes from observations of significant large- 
scale sub-str ucture in the form of stellar streams in the 
halo of M31 lllbata et al.ll2 00l'.'2004t lFerguson et ai]l2n02t 
iMcConnachie et alJ 1200^ .2004) . Other disk-dominated 
galaxies show similar indications of halo sub-structure, 
pointing to a hiera rchical assembly of these systems (e.g. 
IZibetti et allbOO^ . 

Since accretion events tend in general to heat stellar 
disks, a good indicator for a puffed-up disk population is 
the ratio of disk scale length, d, to scale height zq. A minor 
merger with a mass fraction of 0.05-0.2 is expected to ag- 
itate the stellar disk and increase its vertical scale height 
by a factor 1.5-2, compared t o an undisturbed evolving 
disk llVelazauez fc White! ITool . This is confirmed by ob- 
servations o f edge-on spirals with and wi t hout signs of 
interaction llSchwarzkopf fc Dettmarl 120001 1200 ll) . where 
interacting spirals show d/zg < 7. Large-area star counts 
in M31 find that the scale height of the disk varies between 
50 a nd 400 pc, while the scale length is ^ 5—7 kpc ( Hopperl 
Il995i|) . The resulting d/za « 12.5 - 140 is placing the M31 
disk well within the regime of undisturbed stellar disks, 
implying a very small past merging rate. If, however, the 
mass accretion is smooth enough and/or occurs on prefer- 
entially tangential orbits so that the thin-disk integrity is 
unaffected, thin-disk evolution and mass accretion might 
be two spatial ly uncorrelated, but temporally coexisting 
processes fe.g. iFont et allEoOll) . 

It is not clear if some M31 globular cluster s are associ- 
ated w ith t he two stellar st reams discovered bv llbata et alJ 
l|200lh and lMcConnachie et al. ( 2004il . beca use of the still 
too sp arse kinematic information. However. IPerrett et al.l 
l)2003|) find kinematically coupled groups in the M31 glob- 
ular cluster system. Although the location of some globu- 
lar clusters coincides with the general locus of the streams 
(see the northern stream in Fig. it is unlikely that ac- 
cretion contributed significantly to the assembly of the 



M31 globular cluster system, in particular the thin-disk 
sub-population. 

iMorrison et al.l (|22^) find a velocity dispersion in the 
thin-disk globular cluster system of 20 — 40 km/s. 
Globular cluster systems in dwarf galaxies have typical 
velocity dispersions 30 — 80 km/s fe.g. .Olsen et al. 200^ . 
In general, the thin-disk globular clusters in M31 have on 
average too small a velocity dispersion to be accreted from 
merging satellite galaxies, although the very unlikely case 
of a merging trajectory which is perfectly aligned with the 
disk plane could be a way to preserve a small velocity dis- 
persion of the infalling cluster system during the accretion 
event. 

Using the numbers from llbata et al ] ll200lh . the most 
massive stream is estimated to carry 10* to 10^ Mq 
with a total My = —14 mag. Assuming an average 
specific fre quency Sm ^ 10 fo r dwarf galaxies of this 
luminosity i Durrell et alJ Il99fil) . the estimated number 
of accreted globular clus ters is ^ 4. The other stream 
llMcConnachie et al .120041) carries less mass and will there- 
fore contribute even less globular clusters. However, if 
both streams are associated with the two satellites M32 
and NGC 205 (both with a total luminosity My « -16.4 
mag) , the contributed number of clusters might be as large 
as 20-40 objects per galaxy, given typical specific frequen- 
cies Sn ~ 2 — 8 for thes e systems. M32 hosts no globular 
clusters IjAshman fc ZeDf..l998;) . although at least ~ 20 
are expected. 

If the two massive stellar streams added some globu- 
lar clusters to M31, it is unlikely that they are now part 
of M31's thin-disk globular cluster population, based on 
their trajectories. Both stellar streams have very radial 
orbits and therefore do not interact si gnificantly with th e 
thin-disk globular-cluster population jlbata et al.ll200^ . 
Despite their potential past interaction(s), the evolution 
of a kinematically cool thin-disk globular-cluster popula- 
tion could have remained undisturbed. 



5.2. Spatial Distribution 

We find spatially clumped globular cluster sub-popu- 
lations in M31. Intermediate-age globular clusters pref- 
erentially reside inside ~ 5 kpc galactocentric distance, 
which is suggestive of a slightly younger inner-disk/bulge 
cluster population. Clusters with ages below 5 Gyr are 
predominantly found at large radii along the major axis. 
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Most of these globular clusters show thin-disk kinematics 
and have metaUicities below —0.8 dex (see also Fig. [TT|) . 
Their projected positions (see Fig. ^ coincide with lumi- 
nous star- forming regions in M31's outer-disk, that were 
recently observed by GALEX^^. 

This spatial clustering suggests that halo, inner- 
disk/bulge, and thin-disk globular cluster sub-populations 
experienced different formation and/or evolution histo- 
ries. While the intermediate-age population could have 
formed from previously enriched material several Gyr af- 
ter the old halo population was formed, the young clus- 
ter population might be the offspring of star-formation 
triggered by a recent infall of fresh material in the out- 
skirts of M31. However, a direct connection between these 
globular cluster sub-populations and the two giant stel- 
lar streams in the halo is rather unlikely, although we 
find two young metal-rich young globular clusters (B337- 
G68 and B324-S51) close in projection to NGC 205. Since 
there is much more sub-structure visible in M31's halo 
than t he two reported stellar str eams fe.g. lFereuson et aD 
l2002t iMcGonnachie et al.ll20o4 . direct connections be- 
tween other stellar streams and globular cluster sub- 
samples are conceivable. Only a careful kinematic analysis 
of the stellar streams can clarify how much accretion con- 
tributed to the formation of M31's halo and disk globular 
clus ter system. Such a mbitious surveys are on the way (see 
e.g. Ilbata et al.ll2004 . 

5.3. Chemical Enrichment 

In the following we discuss chemical characteristics of M31 
globular clusters that might be informative of their for- 
mation and evolution. While type-II supernovae predom- 
inantly eject a-elements on very short timescales (^ 10 
Myr), the ejecta of type-la supernovae enrich the inter- 
stellar medium with Fe-peak elements with a delay of typ- 
ically - 1 Gvr llGreggiolll997t iMatteucci fc RecchilEoOll 



^ 1- nr^ i II I r ■ -, , ^ 

but see also [Tsui imoto 200jfor a short-living type- la pro- 
genitor). These supernova types mark the two extremes 
of injection timescales for processed material back into 
the interstellar medium. Other elements such as carbon 
and nitrogen are predominantly produced by massive stars 
with a strong depen dence on meta llicity, on timescales 
~ 10-100 Myr (e.g. lMaedejjll99.^. In general, the rel- 
ative delay in progenitor lifetimes between type-II and la 
supernovae implies that abundance ratios can be used as 
tools to clock star-formation timescales. 

A detailed model of the chemo-dynamical evolution of 
globular cluster systems is not available yet and a compar- 
ison of our data with theoretical predictio ns is not possi- 
ble (but see lLi. Mac Low, fc KlessCTll2004 . However, our 
data on [a/Fe] ratios and nitrogen enhancement of M31 
globular clusters allow us to make some qualitative state- 
ments on their formation timescales. 



5.3.1. a-Elements 

Super-solar [a/Fe] ratios of M31 globular clusters can be 
interpreted as the result of relatively short formation peri- 
ods, which ended before a significant number of type-la su- 
pernovae could contribute to chemical enrichment of their 
parent gas clouds. As chemical enrichment increases with 
time it requires successively higher star-formation rates to 
increase the a-element enhancement to a certain level and 
"override" the chemical signature of previous enrichment 
episodes. Hence, the average trend of [a/Fe] is expected 
to be a declining function of metallicity, as enough time 
will be given to type-la supernovae to lower this ratio, un- 
less previous star-formation events were truncated before 
type-la supernovae started to eject large amounts of iron- 
peak elements. In the same line, a decreasing mean [a/Fe] 
is expected at younger stellar-population ages. 

We find no general correlation between [a/Fe] and 
metallicity for our M31 globular cluster sample, which in- 
dicates that the stellar populations in most of these ob- 
jects were predominantly enriched by core-collapse (type- 
II) supernovae. In combination with the rather weak 
anti-correlation of age vs. [a/Fe] and age vs. metallic- 
ity, this suggests that old metal-poor globular clusters 
formed from material which experienced a truncated en- 
richment, dominated by type-II supernova ejecta. Metal- 
rich, intermediate-age and young M31 globular clusters, 
on the other hand, formed from gas that received its met- 
als from type-II and la supernovae, perhaps during the 
simmering star- formation processes in M31's disk. This 
indicates that stellar populations in the M31 disk experi- 
enced significant enrichment by type-la supernovae. The 
dominant fraction of M31 globular clusters, however, was 
enriched by type-II supernovae. 

5.3.2. Nitrogen Enhancement 

Nitrogen is produced in two modes. One speaks of a 
primary and secondary nitrogen component, depending 
on whether the seed elements carbon and oxygen are 
synthesized from helium by the same star (primary) 
or were already in place when the star formed (sec- 
ondary). Sources of primary nitrogen are massive, zero- 
metallicity fast rotating stars that suffer a hypernova 
explosion at the e nd of their v ery short life times (typ- 
ically a few Myr: 
l2nmt iHeger et ahl 



Wooslev fc Weaver 1982: iFrver et alJ 



2ffi 



http:/ /www. galex.caltech.edu/popups/gallery-M31.html 



More massive and faster ro- 
tating hypernova progenitors genera lly inject more ni - 
trogen into the interstell ar medium |Heger et al.|l2000|) . 
However, it is noted by iHeger et al.l that even a non- 
rotating 60 M0 star with a shallow entropy gradient can 
produce strong over-abundances of nitrogen. Secondary 
nitrogen is mainly ejected by strong stellar winds of 
intermediat e-mass stars (^4 — 7 Men) during th eir AGB 
phase fe.g. Ivan den Hoek fc GroenewegeiJll997j) . In gen- 
eral, stars with successively higher mass and higher angu- 
lar momentum produce larger am ounts of both a-elements 
and nitrogen fe.g. Wooslev fc Weaver.. 199,^1. 
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Fig. 13. Correlation between the N/C ratio, computed as the ratio of the NH index and the CH(G4300) index vs. 
age, metallicity, and [a/Fe] ratio. M31 globular clusters are indicated by solid dots. Open circles mark Milky Way 
globular clusters. The hatched region in panel (a) shows the typical locus where Galactic globular clusters reside. A 
linear least-square fit in panel (b) to Galactic globular-cluster data is illustrated by a dashed line. The dashed line in 
panel (c) shows a linear fit to M31 globular clusters. 



This implies that hypernova ejecta have both super- 
solar [a/Fe] ratios and enhanced N/C and N/0 ratios, 
as the production of primary nitrogen is based on the 
consumption of the seed elements carbon and oxygen. 
Moreover, [a/Fe], and N/C and N/0 ratios should be pos- 
itivelv correlated. This is different f or AGB-star yields. 
Ivan den Hoek &: GroeneweeeiJ l)l997j) show that the N/C 
ratio is expected to increase slightly as a function of stel- 
lar mass and metallicity, while N/0 is predicted to have 
a much steeper correlation. Hence, we expect asymptoti- 
cally decreasing N/C and N/0 ratios as a function of in- 
creasing metallicity. We also expect to see a correlation of 
increasing [a/Fe] with increasing N/C and N/0 ratios for 
globular clusters that were mainly enriched by hypernova 
ejecta. 

The importance of primary and secondary nitrogen en- 
richment is currently the subject of a lively debate in the 
literature. The need for a primary plateau-enrichment by 
zero-metallicity massive stars was raised by the detection 
of highly super-solar [N/Fe] ratios in ex tremely metal-poor 
[Fe/H] <, -4) Galactic halo stars (e.g. iNorris et aLll200ll 
l2002HChristlieb et alJ20oll2004^ and metal-poor damped 
Lvman-a svstems fe.g. IPettini et aLl E^02') . 

To study N/C as a function of age, metallicity, 
and [g/Fe], w e match our observations with those of 
iBurstein et alJ l|2004,l and find ten matched M31 globular 
clusters for which NH measurements are available. None 
of the Milky Way globular clusters for whic h NH measure- 
ments were performed by iBurstein et alJ have matched 
age, metallicity, and [a/Fe] ratios derived from Lick in- 
dices in this work. We therefore substitute the latter with 
ages and metalli cities derived from colo r-magnitude dia- 
grams taken from lSalaris &: WeissI l)2002j) and De Angeli et 
al. (2004, private communication). CMD metallicities are 
on the Zinn-West metallicity scale and were corrected for 
an assumed average a-element enhancemen t of [a/Fe] = 
-H0.3 dex, using [Z/H] = [Fe/H]-0.8-[a/Fe] l|Trager et alJ 
l2Q0nt) . No [a/Fe] determinations for the matched Galactic 



globular clusters are available to us at the moment, so 
that we have to postpone a comparison of NH and [a/Fe] 
ratios between Milky Way and M31 globular clusters. 

Th e Galactic and M31 globular cluster samples of 
IBurste in ct al. 2004 cover very different metallicity ranges 
(see panel b in Fig. lU^ . While most M31 globular clus- 
ters fall in the range — 1 ^ [Z/H] ^ 0, most of the Milky 
Way globular clusters have metallicities below —1.0 dex. 
Hence, a direct comparison of NH abundances between 
M31 and Milky Way globular clusters becomes less mean- 
ingful. To bypass this drawback we focus on abundance 
ratios instead. Because of their virtually identical dissoci- 
ation energies l)Tomkin fc Lambertlll98^ . the strength of 
the near-UV NH band and the optical CH band around 
4300 A can be used to derive the nitrogen to carbon (N/C) 
abundance ratio fairly accurately. The strength of the op- 
tical CH-band is measured by the G4300 Lick index which 
is available for all matched M31 globular clusters as well 
as Milky Way globular clusters. We therefore refer to the 
NH/CH(G4300) index ratio in the following as the N/C 
ratio. 

Figure^jshows the N/C ratio as a function of globular 
cluster age, metallicity, and [a/Fe]. We find no clear corre- 
lation of N/C vs. age. Since the age resolution in the cov- 
ered age range (~ 8—13 Gyr) is poor, it would be very in- 
teresting to obtain more NH data for intermediate-age and 
young M31 globular clusters. We note that a correlation 
of increasing nitrogen-enhancement with older age is ob- 
served in starburst galaxies ( Este ban &: Pcim bcrt 1995»). If 
nitrogen-enhancement is mainly produced by hypernovae, 
older globular clusters are expected to show on average 
higher N/C ratios^^. 



For completeness, we note that the Galactic globular clus- 
ter NGC 7078 shows a relati vely high N/C ratio (~ 4). This is 



in contrast to the results of iRauch et al.l i2002l) who measure 
a carbon-enhancement and a nitrogen depletion (implying a 
small N/C ratio) in a sdO cluster star. Further high-resolution 
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We find enhanced N/C ratios in metal- rich M31 glob- 
ular clusters compared to their Galactic counterparts at 
similar metallicities (see panel b in Fig. Almost all 
M31 clusters with [Z/H] ^ —0.8 dex appear to have en- 
hanced N/C abundance ratios, while more metal-poor 
clusters fall on a similar relation as Milky Way globular 
clusters (see dashed line in the middle panel of Fig. [T^ . 
This is clear evidence against an enrichment scenario dom- 
inated by hypernovae, which is expected to imprint higher 
N/C ratios at lower metallicities. Instead, the opposite 
is observed and lends support to the notion that metal- 
rich intermediate-mass AGB stars were given enough time 
(~ 250 Myr) to eject and mix secondary nitrogen into the 
interstellar medium, before the N-enhanced M31 globular 
clusters formed. 

In panel (c) of Figure 1131 we find evidence for a 
correlation of the N/C ratio with [a/Fe], in the sense 
that larger a-enhancements are accompanied by higher 
N/C ratios. A dotted line shows a weighted linear least- 
square fit to the data. By itself this result does not fa- 
vor either of the nitrogen-enrichment mechanisms (hy- 
pernovae or intermediate-mass AGB stars). However, a 
significant hypernova enrichment is expected to produce 
an extreme [a/Fe] signature with values^ 1 dex (e.g. 
IWooslev fc Weaveilll99.4 iFrver et allboOll) . which is ob- 
served in some v ery metal-poor Galactic halo stars (e.g. 
lAoki et aD l2002|) . but not in the integrated light of 
nitrogen-enhanced M31 globular clusters. 

5.4. Chemical Coincidence of Globular Cluster 
Sub-Populations 

It is intriguing that young M31 and LMC globular clusters 
and some counterparts in the M81/Sculptor-group sample 
show very similar CN indices at a given metallicity. Two 
scenarios seem plausible to explain the chemical similar- 
ity: 1) the globular cluster populations experienced similar 
enrichment histories, although they formed in galaxies of 
very different mass and morphology or 2) the young M31 
clusters were stripped from accreted dwarf satellites. 

Most of the young M31 globular clusters reside in the 
outskirts of the M31 disk and show radial velocities con- 
sistent with a kinematically cold thin-disk cluster popu- 
lation. However, it is relatively unlikely that an accretion 
event was so well aligned to leave no kinematical trace of 
a heated globular cluster sub-population. This makes the 
ac cretion scenar i o less likely. 

iForbes et al.l lj2004l) compile data on ages and metal- 
licities of about a dozen Galactic globular clusters (~ 10% 
of the entire system) , which are kinematically associated 
with giant stellar streams and dwarf galaxies in the pro- 
cess of accretion^'^. They find a range of ages and metal- 
licities for accreted Galactic globular clusters, from ~ 6 

measurements for more stars in NGC 7078 are necessary to 
resolve this discrepancy. 

Note that the potentially accreted satellite, which might be 
associated with some of the young M31 globular clusters, does 



to 13 Gyr, and metallicities from [Fe/H] « —2 to solar 
values, which does not allow us to distinguish captured 
from in-situ formed globular clusters using only age and 
metallicity. However, these authors argue that the age- 
metallicity relation for dwarf galaxies (stars and clusters) 
are different from the relation of the remaining Milky Way 
globular clusters. Different age-metallicity relations point 
to different chemical enrichment histories. A comparison 
of the age-metallicity relations of the Milky Way and M31 
globular cluster systems will be the subject of a future 
paper, when more accurate data for the Galactic globular 
cluster system become available. At this point, we note 
that a distinct age-metallicity relation of accreted glob- 
ular clusters does not exclude either of the two former 
scenarios (chemical coincidence and accretion). 

5.5. Formation Scenarios 

There are several formation and assembly scenarios for 
globular cluster systems in spiral galaxies, which can be 
ordered in three different categories. Pre-galactic scenar- 
ios envision globular cluster formation in dense metal- 

J ioor clouds befor e the aggregation of their host galaxy 
Peebles fc Dickelll968j) . In proto-galactic scenarios, glob- 
ular clusters form in loose gas clumps during the coales- 
cence o f the host ga laxy in more or less chaotic merging 
events (jSearle fc Zinnl(l978(l . The third scenario pictures 
an a posteriori assembly where a fraction of the clus- 
ter system is added to the in-situ formed globular clus- 
ter system through acc retion of external globular clusters 
from satellite galaxies ()C6te et al.ll200ol) . It is clear that 
all models overlap in many aspects and all present valid 
formation mechanisms for at least some globular cluster 
sub-populations (i.e. halo, disk, bulge) in M31. 

The current data does not allow a clear-cut answer 
on the importance of each mechanism in the specific case 
of M31, since our data does not sample a representative 
fraction of halo globular clusters. However, our study has 
revealed a globular-cluster age distribution which is en- 
tirely inconsistent with the predictions of the pre-galactic 
and somewhat inconsistent with the proto-galactic sce- 
nario. The formation of intermediate-age globular clusters 
appears only plausible in an extended proto-galactic sce- 
nario where the collapse of the gas clumps occurs over an 
extended period of time, i.e. for several Gyr: a scenario 
closely resembling the hierarchical assembly of galaxies. 
Hence, the pre- and proto-galactic models cannot solely 
describe the assembly of the M31 globular cluster sys- 
tem. Because of their spatial concentration towards the 
inner ~ 5 kpc, an accretion scenario also appears less likely 
as a formation mechanism for the intermediate-age sub- 
population. 

It should be stressed that the formation of individual 
globular clusters might be at stark variance with the sug- 
gested formation pictures. For instance, one might pic- 

not necessarily have to leave a luminous stellar component, and 
hence a traceable stellar stream. 



26 



Puzia et al.: New Light on the Formation and Evolution of M31 and its Globular Cluster System 



ture the formation of intermediate-age globular clusters 
from material which was enriched in the disk and subse- 
quently funneled into the core regions through bar insta- 
bilities. Upscaled version s of this scenario are observed in 
other spiral galaxies fe.g. lJoeee et al.ll2002() . The fact that 
most of the intermediate-age globular clusters show ra- 
dial velocities that are consistent with thin-disk kinemat- 
ics lends support to this formation picture. Alternatively, 
one could im agine a two-phase co llapse, similar to the one 
suggested by iForbes et al.1 l(l997t ) for globular cluster sys- 
tems in early-type galaxies. In this scenario, halo globular 
clusters form first from primordial gas during the proto- 
galactic phase. The first generation of stars enriches the 
remaining gas which later settles in a disk and forms the 
second generation of younger and more metal-rich globu- 
lar clusters. 

A fraction of young M31 globular clusters can be ex- 
plained by the accretion scenario, as we find a few objects 
in close projected vicinity of satellite galaxies and giant 
stellar streams (see Fig. The majority of these young 
relatively metal-rich globular clusters, however, resides in 
the outskirts of M31's disk and forms a kinematically cold 
sub-population. As discussed earlier it is rather unlikely, 
but not impossible, that a well-aligned, pro-grade merging 
event stripped these clusters from an LMC-type satellite 
galaxy. It is also possible that we see evidence for a massive 
star-forming ring (see GALEX press-release images^^) in 
which young star-clusters are formed from recently com- 
pressed molecular clouds, which might be triggered by re- 
cent tidal stress. 

For completeness we note that the old globular cluster 
sub-population has a wide range of metallicities and its 
spatial distribution is consistent with a halo population. 

Chemical compositions of intermediate-age and young 
M31 globular clusters are consistent with a formation on 
extended timescales. The chemical signatures of their old 
counterparts are expected to be the result of short star- 
formation events in which the interstellar medium was sig- 
nificantly polluted by the thermonuclear ashes of type-II 
supernovae. In particular, iLi & BursteinI (|2D03l) put for- 
ward the idea that the old M31 globular-cluster pop- 
ulation experienced significant enrichment by nitrogen- 
rich hypernova ejecta, motivated by the assumption that 
the nitrogen overabundance in metal-poor Galactic halo 
stars and old M31 globular clusters has the same origin. 
Although conceptually plausible, we do not subscribe to 
this idea in its details. It is likely that hypernovae con- 
tributed significantly to the enrichment of the interstellar 
medium from which the most metal-poor Galactic stars 
([Z/H] ^ —4) have formed. However, their metallicity is 
about two orders of magnitude lower than that of the most 
metal-poor known Galactic and M31 globular cluster^^. It 
is hard to believe that no moderate-mass 10 Mq) type- 
II supernovae contributed to the chemical signatures of the 



e.g. pittp://antwrp. gsfc.nas a.gov/apod/ap031222. htrnl] 

The most metal-poor Galactic globular cluster NGC 5053 
has a [Fe/H] « -2.3 dex lHarris..l99&'l . 



most metal-poor globular clusters since the early epochs 
of chemical enrichment. 

Furthermo r e, th e hypernova scenario advocated by 
iLi fc Bursteh]l|2003l) . has difficulties in explaining the for- 
mation of the old metal-poor, but not nitrogen-enhanced 
M31 globular clusters (such as B311, see Fig.^| this glob- 
ular cluster has the smallest N/C ratio in our sample). It is 
conceivable that this might be a result of chemical variance 
in the very early epochs of chemical enrichment. Clearly, 
the hypernova scenario requires some detailed modeling 
and comparison with observations. 

It would also be interesting to check if there is a range 
in nitrogen enhancement for Milky Way and other glob- 
ular cluster systems at low metallicities, resembling the 
chemical variance during the very early star-formation 
epochs. We have embarked on a photometric and spectro- 
scopic multi-wavelength study of globular clusters in the 
Milky Way and Local Group galaxies with ground and 
space-based telescopes and will address these issues in the 
future in more detail. 



6. Summary and Conclusions 

The presence of a populous M31 globular-cluster subsys- 
tem significantly younger than found in the Galactic glob- 
ular cluster system is certainly the most surprising result 
of this study. Evidence for the counterpart globular cluster 
population in the Milky Way is missing, but we have found 
for the presence of similar globular cluster sub-populations 
in other spiral galaxies. Previous studies have shown that 
such intermediate-age and young globular clusters also ex- 
ist in giant ear ly-type llPuzia et a l."2002 bl: iPuzis 
iLarsen et a,].l2h n3: Bri dges et a,]^ ,2nn3: 1 
and dwarf galaxies (e.g. iBeaslev et al.l li 

In general, the M31 globular cluster system hosts a 
wider variety of ages and chemical compositions than the 
Galactic globular cluster system, and is therefore more 
representative for extragalactic stellar populations, espe- 
cially for the comparison with stellar populations in early- 
type galaxies. It seems that globular cluster system in 
spiral galaxies might have experienced conceptually sim- 
ilar, but quantitatively different, formation and/or as- 
sembly histories. Some extreme cases might even resem- 
ble globular-cluster formation/assembly histories of early- 
type galaxies. 

Below we summarize the main results of this work. 

• In addition to a dominant population of globular clus- 
ters with ages > 10 Gyr and metallicities from [Z/H] « 
—2.0 to solar values, we find a population (~ 20 ± 7% 
of the observed sample) of intermediate-age globular 
clusters with ages between ~ 5 and 9 Gyr and a mean 
metallicity around —0.6 dex. A few of these objects 
have metallicities around —1.6 dex. 

• We confirm the presence of young M31 globul ar clus- 
ters recently identified bvlBeas lev et alJ ((2004 ). which 
have ages ^ 1 Gyr and relatively high metallicities 
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around —0.4 dex; however, two metal-poor young ob- 
jects have also been found. 

• The M31 globular cluster system has a clearly super- 
solar mean a/Fe abundance ratio of -1-0.14 ± 0.04 dex 
with a dispersion a = 0.32 dex. However, the distri- 
bution appears bimodal with peaks around -1-0.2 and 
—0.03 dex, which have different dispersions, 0.19 and 
0.50 dex, respectively. 

• Intermediate-age and young globular clusters show 
roughly solar abundance ratios. 

• We find evidence for an age- [a/Fe] relation in the sense 
that younger clusters have smaller mean [a/Fe] ratios. 
This is driven by the metal-rich intermediate-age and 
young cluster population with metallicities higher than 
[Z/H] « -1.0 dex. 

• The indices of globular clusters older than ^ 5 Gyr are 
consistent with a factor of three or higher in nitrogen- 
abundance enhancement compared to their younger 
counterparts. 

• Using kinematical data for a thin-disk model 
I Morrison et al.l Eo04|l we find that the halo popula- 
tion of globular clusters is old 9 Gyr), has a bi- 
modal metallicity distribution and super solar a/Fe- 
enhancement. 

• Thin/ thick-disk globular clusters have a wide range of 
ages and are predominantly metal-rich with a slightly 
smaller mean [a/Fe] ratio. 

• With the structu ral parameters of M3 1 globular clus- 
ters measured bv iBarmbv et al.l iooi), we confirm a 
previously found correlation between half-light/tidal 
radius and globular clusters metallicity, which is most 
likely due to the correlation of half-light/tidal radius 
and galactocentric distance. 

• Comparing the integrated light of the M31 and Milky 
Way bulge, wc find that both host metal-rich, a- 
element enhanced and old stellar populations. 

• We find indications for an enhancement in carbon and 
nitrogen in the central region of M31, while the inte- 
grated light of outer M31 stellar populations is consis- 
tent with a nitrogen enhancement only. 
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T able A.l. Lick i ndices for our sample M31 globular clusters. Indices were comp uted with the passband defin itions 
of IWorthev et all l)l994|) an d for the higher-ord er Balmer lines the definitions of IWorthev fc: Ottavianil l)l997f) . We 
adopted the nomenclature of iGalleti et alJ l|2004j) for all our objects. Note that B166 is suspected by some studies to 
be a foreground star. Since the case is not clear, we mark this object as potentially misclassified but keep it in our 
globular cluster sample. However, we point out that all other objects we re classified by seve ral independent studies as 
genuine globular clusters, based on their colors and radial velocities (see iGalleti et al]l2n04l) . 
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Table A. 3. Ages, Metallicities, and [a/Fe] ratios for M31 globular clusters. We point out that the quoted uncertainties 
are the propagated statistical uncertainties of the data in combination with the SSP models used in this work. 
Systematic uncertainties of the presented parameters are of the order (for [Z/H] and [a/Fe]) or larger (for age) than 
the statistical errors and are discussed in Section 3.3. 
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